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Abstract 

The  method  of  matched  asymptotic  expansions 
is  used  to  analyze  the  interaction  of  a normal 
shock  wave  with  an  unseparated  turbulent  bound- 
ary layer  on  a flat  surface  at  transonic  speeds. 

The  theory  leads  to  a three-layer  description  of 
the  interaction  in  the  double  limit  of  Reynolds 
somber  approaching  infinity  and  Mach  number  ap- 
proaching unity.  The  interaction  involves  an  out- 
er, Imrlscid  rotational  layer,  a constant  shear- 
stress  wall  layer,  and  a binding  region  between 
them.^The  pressure  dist  utlcn  Is  obtained  from 
a numerical  solution  of  th.  outer-layer  equations 
by  a mixed-flow  relaxation  procedure.  An  analyt- 
ic solution  for  the  skin  friction  is  determined 
from  the  inner-layer  equations.  The  significance 
of  the  mathematical  model  is  discussed  with  ref- 
erence to  existing  experimental  data. 

Nomenclature 

a = /m"#/l  + m# 

at  constant  in  Bradshaw’s  model,  usually 

0.15 

B Law  of  the  wall  profile  parameter 
C,  skin  friction  coefficient,  - t * /\  ptu£* 

C,  function  determined  from  blending-layer 

solution,  eq.  (5,26) 

D damping  factor,  eq.  (6.1) 

c turbulcrt  kinetic  cnrigy 

E perturbation  turbulent  kinetic  energy  in 

blending  layer 

F wall  layer  Reynold,-,  stress  function,  eq. 

(6. 13) 

g,  defined  try  eq.  (A.  7) 

H defined  by  cq.  (C>.  14) 

K viscosity  function,  cq.  1G.9) 

„ M'?  - 1 . 

K|  = ' Vm*~'  ’ 'r*in‘s,orni,-f'  interac t>on  pa  ram - 
* * otor 

Kt  ^K,/>  4 1 , interaction  similarity  parameter 
L reference  length 

dissipation  length 
M Mach  number 


p pressure 


p - y—-  Pi  t transformed  perturbation  pres- 
Xt  sure  In  outer  region 
P perturbation  pressure  in  blending  layer 


R wall  layer  density  function,  eq.  (6. 14) 

ReL  ■ pi  u*  L*/pi , Reynolds  number 
Re^  «pi*u*6|/pi,  Reynolds  number  based  on 
boundary-layer  thickness  upstream  of 
the  interaction 
T temperature 

u velocity  in  the  X direction 

u*  wall  layer  velocity  function,  eq.  (6.4) 
ut  defect* velocity  In  the  outer  region  up- 
stream of  the  interaction 
U|  u*  evaluated  upstream  of  the  interaction 

fi  Van  Driest  generalized  velocity,  eq.  (3.1), 


U perturbation  velocity  in  x direction  in  the 
blending  layer 

uT  =■  'TYJpy,  , friction  velocity 

v velocity  in  the  y direction 

V perturbation  velocity  in  y direction  in  the 

blending  layer 

x*  coordinate  parallel  to  the  plate 

x - x*/L* 

x stretch  x coordinate  defined  in  eq.  (4.1) 

* */^xI 

y*  coordinate  normal  to  the  plate 

? r y*/L* 

y stretch  y coordinate  in  outer  region,  eq. 

(4.2) 

Y stretched  y coordinate  in  hlumling  layer, 

eq.  (5.G) 

y’  stretched  V coordinate  in  the  Inner  layer, 
cq.  (3.9) 

a Bradshaw’s  mean  dilatation  parameter 

y ratio  of  specific  heats,  ~ i.  4 


^Thc  latter  plume  of  this  research  was  partially  supported  by  the  Office  of  Naval  Research  under  Contract 
N00014-71-C4636.  i’ome  of  the  ideas  regarding  the  turbulence  structure  were  developed  in  conjunction 
with  NASA  Contract  NAS  1-12426. 

•Director  of  Acrophyslcs,  Research  Department 
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r perturbation  Reynolds  stress  In  Mending 

layer 

6 boundary  layer  thickness 

A,  = Ml  - 1,  shock  strength  parameter 

A,  x length  scale 

A,  y length  scale 

« * JCt  J2,  nondimensional  friction  velocity 

upstream  of  the  interaction 
| » l/c*ReL 

t4  turbulent  dissipation  function 

rj  y coordinate  in  the  computational  plane, 

eq.  (7.16) 

k Von  Karman  constant,  - 0.  41 

X *-(x»  f2ut+p,)/y+l 

* = VXt 

H coefficient  of  viscosity 

$ x coordinate  in  the  computational  plane, 

eq.  (7.15) 

t law  of  the  wake  profile  parameter 

p density 

p,  = 2H1.U, , defect  density  in  outer  layer  up- 

stream of  the  Interaction 
o defined  in  eq.  (5.23) 

r Reynolds  stress 

<J>  velocity  potential  for  first-order  outer 

solution,  eqs.  (4. 12),  (4. 13) 
d>tft  split  velocity  potential,  eqs.  (7.11),  (7.12) 

<5  =(y+l)<5/Xt/l 

Xt  = Mi  - 1/c  Interaction  parameter 

u relaxation  parameter 

(i  generalized  potential  function  in  the  blend- 

ing layer,  cq.  (5.29),  eq.  (5.30) 

Subscripts 

« free  stream  conditions 

0 conditions  upstream  of  the  interaction/ 

w wall  condition  / 

1, 2 perturbation  indices  j 

Superscripts 

* dimensional  quantity 

* turbulent  fluctuation 

1.  Introduction 

Interactions  between  shock  waves  and  turbu- 
lent boundary  layers  occur  in  a number  of  practi- 
c'd aerodynamic  situations.  Examples  include 
supercritical  flows  over  airfoils,  exhaust  nozzles, 
and  inlets.  These  interactions  often  lead  to  pre- 

mature separation  that  can  cause  severe  flow 
problems  such  as  buffeting,  inlet  flow-field  dis- 
tortion, and  excessive  boattail  drag.  Shock  wave- 
boundary  layer  Interactions  are  also  known  to  be 
a major  source  of  significant  scale  (Reynolds 

number)  effects  observed  in  transonic  wind-tun- 

nel tests. 


In  the  present  investigation,  we  apply  asymp- 
totic and  numerical  methods  to  the  solution  of  the 
interaction  of  a normal  shock  wave  with  a turbu- 
lent boundary  layer.  We  restrict  our  analysis  to 
the  simplest  case  of  a weak  shock  wave  impinging 
upon  a smooth  flat  plate.  The  shock  strength  is 
assumed  to  be  -sufficiently  sma’l  as  to  preclude 
boundary- layer  separation,  and  for  standard  tran- 
sonic approximations  to  apply  in  the  invisetd  part 
of  the  flow  field.  We  also  assume  that  the  Rey- 
nolds number  Is  large  apd  that  the  profile  of  the 
approaching  turbulent  Ixjundary  layer  is  fully  de- 
veloped. Homme  of  these  assumptions,  the  pres- 
ent study  should  bo  viewed  as  a Hist  step  towards 
the  solution  of  more  important  problems  involving 
separation,  stronger  shock  waves,  and  flows  over 
curved  surfaces.  It  is  hoped  that  the  present  ef- 
fort will  lead  to  a iheoretical  framework  that  can 
eventually  include  these  effects. 

Although  numerous  experimental  and  analytical 
investigations  of  turbulent  interactions  lave  been 
conducted  over  the  past  few  decades,  our  under- 
standing ot  the  general  problem  area  remains  in- 
complete, even  for  the  basic  case  considered  here 
(see  Ref.  1 for  a review  of  the  shockwave-turbu- 
lent boundary  layer  interaction  literature). 

A basis  for  a rational  treatment  of  turbulent 
interacting  Hows  was  established  in  a pioneering 
study  by  L.iKhthilll * * **>  in  1953.  In  that  work,  a two- 
layer  model  of  an  Interacting  turbulent  boundary 
layer  was  developed,  which  was  based  on  preced- 
ing efforts  by  Howarth/5*  and  Tsicn  aqd  Finston/*1 
Lighthill’s  model  assumed  that  the  streamwise 
length  scale  Is  small  at  high  Reynolds  numbers, 
so  that  molecular  and  turbulent  stresses  could  be 
neglected  over  most  of  the  boundary  layer.  This 
leads  to  a two-layer  description  of  an  interacting 
turbulent  boundary  layer:  an  outer  inviscid,  rota- 
tional layer  where  normal  pressure  gradients  are 
important,  and  an  inner  viscous  layer  required  to 
satisfy  the  no-slip  condition  at  the  surface.  Al- 
though the  model  was  proposed  by  Lighthill  over 
20  years  ago,  the  structure  of  the  inner  layer  and 
details  of  the  matching  of  the  inner  and  outer  solu- 
tions have  not  been  clarified.  It  is  interesting  to 
note  that  the  general  ideas  of  the  Lighthill  model 
served  as  a basis  for  the  development  of  rational, 
asymptotic  methods  for  laminar  flows  (e.g.,  Refs. 
5-10). 

Recently,  the  two-layer  model  was  employed 
by  Roshko  and  Thomkc,m>  and  Elfstromtu>  to 
treat  the  supersonic  turbulent  interaction  In  a 
wedge  compression  corner.  About  the  sumo  time, 
Watson,  Murphy,  and Rose,llJ> and Rosc<,w applied 
the  model  to  Interaction  with  oblique  shock  waves. 
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In  Ref.  11  and  13  calculations,  the  Inner  layer  was 
Ignored  and  the  solution  of  the  outer,  rotational 
inviscid  byer  was  solved  by  the  r.'ethod  of  charac- 
teristics. (Reference  12  used  a semi- analytic 
method  to  solve  the  supersonic  region.)  In  all 
these  calculations,  the  Initial  profile  had  to  be  cut 
off  arbitrarily  to  ensure  that  the  flow  remained 
supersonic  in  the  outer  region.  Experimental  data 
presented  in  these  studies  clearly  Indicate  that  the 
rise  in  surface  pressure  occurs  in  two  stages. 
First,  the  pressure  rises  nearly  discontlnuously 
to  a clearly  defined  level  over  a streamwlse  dis- 
tance of  a boundary- layer  thickness  or  less.  This 
is  followed  by  a gradual  increase  in  pressure  to 
an  asymptotic  value  far  downstream  of  the  inter- 
action. In  some  cases  overshoots  were  observed. 
The  computations  treated  the  initial  pressure  rise 
as  a discontinuity.  Comparisons  with  the  experi- 
ments Indicated  that  good  agreement  for  the  sur- 
face-pressure distribution  could  be  obtained  If  the 
cutoff  was  chosen  to  match  the  pressure  at  the 
downstream  end  of  the  initial  pressure  rise. 

These  results  give  support  to  the  main  features  of 
the  two- layer  model  and  demonstrates  the  useful- 
ness of  representing  the  outer  part  of  the  boundary 
byer  by  a rotational  Inviscid  flow.  However,  the 
two-layer  model  in  its  present  form  Is  incomplete. 
It  does  not  treat  the  flow  in  the  initial  pressure 
rise,  nor  consider  the  structure  of  the  flow  in  the 
inner  layer  and  its  elfect  on  the  outer  solution. 

As  a result  there  is  little  indication  of  how  the 
slip  velocity  should,  even  In  principle,  be  deter- 
mined. 

For  the  transonic  Interactions  being  considered 
In  the  present  study,  the  overall  pressure  rise  is 
more  gradual  with  the  result  that  it  ran  be  de- 
scrilicd  using  a single  st Teamwise  length  scale. 
This  is,  essentially,  because  the  sonic  line  occurs 
In  the  main  part  of  the  boundary  layer,  well  away 
from  the  wait,  when  the  free  stream  Mach  number 
is  close  to  one.  We  will  lx*  able  to  oUaln  a com- 
plete solution  lor  the  surface  pressure  distribu- 
tion without  an  c mpiricu)  velocity  cutoff. 

The  presold  analysis  in  based  on  the  application 
of  the  method  of  matched  asymptotic  expansions 
to  the  full  Na vice- Stokes  equations  in  the  limit  of 
large  Reynolds  number  and  Mach  number  ap- 
proaching one  (weak,  normal  shock  waves).  We 
show  that  these  limits  lead  to  a three-layer  de- 
scription of  th'*  boundary  layer  including  an  invis- 
cid outer  region,  a dissipative  inner  region,  and 
a blending  region  iietwecn  them.  The  analysis 
can  bo  viewed  as  a natural  extension  of  the  asymp- 
totic theories  of  Mcllor,mi  Yajiiik,m>  Bush  and 
Fendcnm,u*  for  incompressible,  and  Afzal<m  for 
compressible  non- interacting  turbulent  Iwundary 


layers.  These  authors  show  that  the  conventional 
defect  and  bw  of  the  wall  description  of  a turbu- 
lent profile  arc  related  to  asymptotic  solutions  of 
the  Navler-Stokes  equations. 

Our  theoretical  model  differs  from  Light  hill’s 
In  the  defect  (outer)  region  mainly  due  to  appear- 
ance of  an  additional  non-linear  term  in  the  equa- 
tion of  motion  due  to  transonic  considerations.  In 
addition,  we  note  that  the  inltbl  profile  can  be  ex- 
panded in  small  defect  form  for  brge  Reynolds 
numbers.  This  significantly  simplifies  the  prob- 
lem, since  it  allows  us  to  represent  the  solution 
In  the  defect  region  as  a small  perturbation  to  a 
weakly  sheared  oncoming  stream.  A velocity  po- 
tential can  be  introduced  and  the  governing  equa- 
tion can  be  written  in  a form  suitable  for  numeri- 
cal solution  by  mixed  flow  relaxation  procedures, 
first  developed  by  Murman  and  Cole.<M> 

A similar  approach  to  turbulent  interaction 
problems  based  on  tl  s method  of  matched  asymp- 
totic expansions  has  been  carried  out  by  Adamson 
and  Fco!l,,  Their  results  apply  to  much  weaker 
shock  waves  than  those  consideied  in  the  present 
study.  They  applied  the  theory  to  interactions  of 
very  weak  oblique  shock  waves  and  sho  wed  that 
normal  pressure  gradients  could  be  ignored  in  the 
defect  byer.  (This  is  not  true  for  the  stronger 
shocks  considered  here.)  They  also  derived  an 
analytic  expression  for  the  solution  to  a turbulent 
free  interaction.  In  our  presentation,  we  show 
that  our  outer-layer  equations  reduce  to  those 
of  Ref.  21  in  the  limit  as  a normalized  shock 
strength  parameter  approaches  zero. The  descrip- 
tion of  tlx*  inner  layer  differs  only  in  that  Adam- 
son and  Feo  employ  an  eddy  viscosity  model  of  the 
Reynolds  stress,  whereas  the  present  work  uses 
a turbulent  kinetic  energy  approach.  The  results 
of  our  analysis  agree  with  their  estimate  of  the 
pressure  rise  for  incipient  separation  and  that  thi: 
pressure  rise  should  uc  independent  of  Reynolds 
numlier  (lor  tin*  large  Reynolds  numbers  con- 
sidered). 

We  note  that  lnger  and  Mason<m  have  recently 
considered  the  extension  of  LighthlU’s  model  to 
transonic  flows  and  have  called  attention  to  (he 
importance  of  Including  the  non-linear  transonic 
term  In  the  equation  of  motion. 

In  the  following  section,  wo  set  out  the  govern- 
ing equations  and  discuss  the  considerations  used 
to  establish  the  order  of  magnitude  of  the  mean 
flow  variables  and  turbulent  correlations.  The 
structure  of  the  solution  in  the  double  limit  of 
Mach  numlicr  going  to  one  and  Reynolds  number 
approaching  infinity,  along  with  the  length  scales 


appropriate  to  the  various  regions  will  be  estab- 
lished. Section  HI  considers  the  solution  in  the 
non-interacting  part  of  the  boundary  layer  up- 
stream of  the  shock  wave.  Here  we  will  show  that 
the  results  obtained  from  the  method  of  matched 
asymptotic  expansions  for  incompressible  flow 
must  be  considerably  modified  to  account  for  com- 
pressibility. Section  IV  considers  the  outer  or 
defect  layer.  It  will  be  shown  that  the  flow  in  this 
region  is  governed  by  a generalization  of  the  usual 
non-linear  transonic  potential  flow  equation  and 
that  the  solution  depends  on  two  similarity  param- 
eters. Next,  the  basic  equations  governing  the 
blending  layer  are  derived  in  Section  V.  We  show 
that  history  effects  arc  Important  in  this  region 
and  suggest  specific  turbulence  models  foi  closing 
the  system  of  equations.  In  Section  VI,  we  de- 
velop the  solution  in  the  wall  layer  and  carry  out 
the  matching  with  the  blending  layer,  A solution 
for  the  skin  friction  is  then  obtained.  The  com- 
plete boundary  value  problem  for  I ho  numerical 
solution  of  the  outer  region  is  formulated  in  Sec- 
tion VII.  The  numerical  methods  used  to  solve 
the  problem  are  discussed  in  detail.  Section  VIII 
contains  results  of  the  numerical  solution  fo.' 
several  normal  shock-wave  interaction  problems 
along  with  comparisons  and  discussions  of  experi- 
mental data.  Finally,  a summary  and  a critique 
of  the  major  findings  of  this  effort,  along  with  a 
discussion  of  methods  for  improving  and  extending 
the  analysis,  are  presented  in  Section  IX. 


II.  Formulation 

The  flow  field  under  consideration  is  the  inter- 
action of  a weak  normal  shock  wave  with  a turbu- 
lent boundary  layer  on  a flat  surface  as  sketched 
in  Fig.  1.  For  simplicity,  the  wall  is  considered 
to  be  smooth,  adiabatic,  and  non-porous.  We  as- 
sume the  interaction  to  occur  at  a distance  L* 
from  the  leading  edge  of  the  plate  and  that  the  flow 
is  uniform  upstream  of  the  interaction,  outside 
the  boundary  layer  with  velocity  ut,  density  p!, 
pressure  />!  and  viscosity  pt.  (The  asterisk  de- 
notes dimensional  quantities.)  These  quantities 
define  a Reynolds  number 


L n; 


The  shock  wave  is  assumed  to  lx;  normal  to  the 
flat  plate  far  above  the  Interaction  and  can  bo 
characterized  by  the  value  of  a single  shock 
strength  parameter,  which  we  take  to  be 


where  M„  is  the  uniform  Mach  number  upstream 
of  the  shock. 

Limit  Process 

We  will  consider  the  solution  of  this  problem  in 
the  double  limit  of 

Reb-«>  and  <1,-0  (2.3) 

It  lias  been  found  to  be  more  convenient  to  arrange 
the  expansion  in  terms  of parameter  i,  equal  to 
the  non-dimeusinual  friction  velocity  of  the  ap- 
proaching boundary  layer,  defined  as 

< » /Cf/2  (2.4) 

where  C.  is  the  skin  friction  coefficient.  Since  € 
may  bo  related  to  the  Reynolds  number  by  the  so  - 
lution ol  the  undisturbed  boundary  layer,  tlie  limit 
defined  in  eq.  (2.  3)  may  be  equivalently  written  as 

€ — 0 and  A,-0  (2.5) 

In  the  limits  given  in  eq.  (2.5)  (or  (2.3)),  the 
boundary- layer  thickness  vanishes  and  the  shock 
wave  becomes  infinitesimally  weak.  The  standard 
experimental  correlations  of  low-speed  turbulent 
boundary-layer  data  (e.g.,  see  Ref.  23),  and  the 
more  recent  asymptotic  analyses, show  that 
the  velocity  profile  lias  a small  defect  form  at  high 
Reynolds  numbers.  The  velocity  differs  from  free 
stream  by  a term  of  0(c)  over  most  of  the  bound- 
ary layer.  The  small  defect  form  of  the  turbulent 
profile  is  an  important  element  in  our  approach. 
We  will  subsequently  show  that  the  two  limits  in 
eq.  (2. 5)  are  not  independent  in  tlie  most  general 
case  and  that  a distinguished  limit  does  exist. 

This  limit  can  best  be  defined  in  terms  of  a simi- 
larity variable,  x%  defined  by 


The  limit  is  then  assumed  to  be 


f-0  } 

M.-l}*1  "* 


In  this  limit,  the  velocity  jump  across  the 
shoc  k wave  is  the  same  order  as  the  velocity  de- 
tect in  the  approaching  boundary  layer.  In  the  in- 
teraction problem  of  Adamson  and  Foo,<iu  y,  van- 
ishes and  the  shock  strength  is  weaker  than  con- 
sidered here.  At  the  other  extreme  the  present 
results,  ami  those  in  Ref.  21,  indicate  that  the  shock 
strength  for  inripent  separation  is  0(1),  indepen- 
dent of  Reynolds  number.  Thus,  it  follows  that 
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for  separation  x,  - Oic*1)-  « which  is  stronger  than 
the  shock  strength  considered  here. 

Length  scales 

For  the  sake  oi  clarity  *n  this  section,  we  will 
describe  the  asymptotic  structure  oi  the  flow  field 
aod  will  give  the  pertlne  it  length  scales  in  terms 
of  < without  mathematics,  justification.  The  proof 
that  the  suggested  length  scales  are  appropriate 
naturally  follows  from  the  substitution  of  the  as- 
sumed expansions  into  the  equation  of  motion  and 
from  a careful  comparison  of  the  orders  of  magni- 
tude of  the  various  terms.  This  part  of  the  analy- 
sis Is  conventional  and  we  will  not  comment  fur- 
ther upon  it. 

At  large  Reynolds  numbers,  the  boundary  lay- 
er spreads  the  discontinuous  pressure  rise  Im- 
posed by  the  impinging  shock  wave  over  a 
small  streamwise  distance.  There  has  been  some 
uncertainty  In  previous  studies  of  turbulent  inter- 
actions about  the  magnitude  of  this  interaction 
length.  In  the  limit  situation  defined  by  eq.  (2.  7), 
it  can  be  shown  that  the  flow  lie  haves  as  an  invis- 
ctd  fluid  over  most  of  the  boundary  layer  and  that 
the  sonic  line  lies  within  the  main  portion  of  the 
boundary  layer.  The  vertical  exl.na  of  the  Inter- 
action zone  is  fixed  by  the  upstream  flow  to  be*  of 
the  order  of  a boundary  layer  thickness.  For  a 
fully  developed  turbulent  boundary  layer,  the 
thickness  scales  with  the  frkiloti  velocity,  and 
thus 

5J/L*  - OU)  (2.6) 

where  Is  the  boundary  layer  thickness  up- 
stream of  the  interaction.  With  velocity  distur- 
bances of  0(c),  it  follows  front  conventional  tran- 
sonic con  side  rations  that  the  streamwise  length 
scale  Is  given  by 

v 0(*m)L*  (2.0a) 

As  men* lotted  above,  the  y ordinate  scales  as 

(Aj)ouUr  6S  0(c)L*  (2.0b) 

Note  that  the  order  of  the  streamwise  length  scale 
is  less  than  a loundary  layer  thickness.  Physical- 
ly, this  follows  from  the  observation  that  the  ex- 
tent of  upstream  influence  is  controlled  by  the  lo- 
cation of  the  sonic  line  in  the  boundary  layer.  For 
the  present  case,  consideration  of  the  slo|>c  of  the 
characters  sties  in  the  upstream  flow  leads  direct- 
ly to  eq.  (2.0). 


The  length  scales  defined  in  eq.  (2. 9)  describe 
the  outer  or  defect  region  of  the  interaction, 
shown  schematically  In  Fig.  2.  Under  the  limit 
In  eq.  (2. 7)  and  the  above  length  scales,  all  turbu- 
lent and  viscous  stresses  are  negligible.  Hence, 
as  In  Llghthill’s  ad  hoc  model,  tlie  flow  in  this 
region  is  governed  try  rotational  invisetd  flow 
equations. 

The  solution  of  the  invlscld  outer-layer  equa- 
tions does  not  satisfy  the  no-slip  condition  at  the 
surface  and,  hence,  at  least  one  additional  inner 
region  is  required.  The  present  analysis  shows 
that  two  inner  regions  are  required  (as  shown  in 
Fig.  2).  One  is  a conventional  compressible  wall 
layer.  In  this  region  the  flew  is  in  local  equilib- 
rium, where  the  solution  for  the  shear  stress  ad- 
justs instantaneously,  in  response  to  changes  in 
velocity  at  the  outer  edge.  Thus,  the  impinging 
shock  wave  Interacts  with  the  main  part  of  the 
boundary  layer  and  decelerates  ’.he  flow  at  the 
outer  edge  of  the  wall  layer.  This  reduction  in 
velocity  causes  a reduction  in  t tear  stress  and 
skin  friction  in  the  wall  layer.  *>om  a considera- 
tion of  the  turbulent  klnetlc-energ;  equation,  it 
can  be  demonstrated  that  the  shear  stress  in  the 
cuter  part  of  the  boundary  layer  is  “frozen"  at 
the  upstream  values  (to  lowest  order  in  *).  Hence, 
a discontinuity  develops  between  the  shear  stress 
in  the  outer  aiul  wail  layers.  The  resolution  of 
this  mismatch  requires  the  introduction  of  a third 
region  called  the  blending  layer.  This  three-lay- 
er structure  appears  to  be  characteristic  of  tur- 
bulor.t  boundary  layers  in  a steep  pressure  gradi- 
ent. 

The  length  scale  for  the  wall  layer  is  *»ven  in 
terms  of  the  local  friction  velocity  in  t\e  usual 
manner  by 

(AvUl  l»r*r : p*(x)/p*(x)u*(x)-0(<c)  (2. 10) 

where  pj  is  the  local  viscosity  coefficient,  pj(x) 
the  local  density,  and  u*the  local  friction  velocity, 
all  defined  with  respect  to  surface  values,  and  c 
is  defined  (following  Melloru$*)  by 

i-'^—0{p*|.(-1A)|  (2.11) 

The  length  scale  in  the  blending  layer  is  deter- 
mined by  introducing  cq.  (2.9a)  into  the  streamwise 
momentum  equation  together  with  the  requirement 
that  inertia  terms  and  turbulent  sijear  stress 
t°rms  Ik?  of  the  same  order  of  magnitude  near  the 
wall.  Tlds  leads  to  the  length  scale  in  the  blend- 
ing layer  t/ 

^ 0(c*'*)L*  (2. 12) 
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Equations  (2.8)  through  (2.12)  provide  the  length 
scales  that  are  to  be  used  to  constrict  asymptotic 
expansions  of  the  solution  in  each  of  the  three  re- 
gions. It  can  be  noted  here  that  to  lowest  order  in 
<,  normal  pressure  gradients  are  negligible  in  the 
wall  and  blending  layers.  Hence,  the  pressure 
distribution  through  the  Interaction  is  obtained 
from  a solution  of  the  outer-layer  equations. 

Also,  it  will  be  shown  that  a solution  for  the  skin 
friction,  to  second  order,  can  be  obtained  from 
the  wall-layer  equations  with  the  outer-layer  sur- 
face pressure  imposed.  These  solutions  will  be 
discussed  In  Sections  IV  and  VI. 

Governing  Equations 

The  present  formulation  is  bused  on  the  time- 
averaged  Navier-Stokes  equations.  Since  turbu- 
lent non-equilibrium  or  history  effects  arc  likely 
important  in  steep  pressure  gradients,  we  will 
use  Bradshaw’s  turbulent  kinetic-energy  formula- 
tlon(*<)  in  our  closure  model.  However,  our 
analysis  will  show  that  the  need  for  a specific 
closure  model  first  arises  in  the  solution  for  the 
second-order  problem  in  the  blending  layer  and 
affects  oidy  the  second-order  solution  for  the  skin 
friction. 

A large  number  of  turbulent  correlations  arise 
In  the  time-averaged  comoresslblc  Navler-Stokes 
equations.  We  will  follow  Bradshaw VI4>  and 
Afzars<,#*  use  of  Kistler’s,|5>  data  to  estimate  the 
order  of  magnitude  of  the  correlations  in  terms 
of  the  undisturbed  friction  velocity  c and  neglect 
the  terms  that  are  of  higher  order  in  all  three  re- 
gimes of  the  present  problem.  In  these  papers, 
it  was  shown  tha*  the  fluctuations  of  velocity, 
te.  >peraturc,  a’ d density  are  of  0(c)  and  :hc  fluc- 
tua  'ons  in  pressure  are  of  0(c*).  In  using  these 
estimates  to  ellm*  vate  certain  higher-order  cor- 
relation terms  we  must,  of  course,  take  note  of 
the  small  streamwise  length  scale  of  the  interac- 
tion. We  also  note  that  the  viscous  terms  are 
only  Important  in  the  innermost  wall  layer  where 
boundary-layer  approximations  arc  valid.  Hence 
only  the  boundary-layer  form  of  the  viscous  terms 
in  the  streamwise  momentum  equation  will  be  ex- 
plicitly included.  Thus,  the  governing  equations 
(from  Ref.  24)  containing  only  those  terms  im- 
portant In  any  of  the  three  layers  arc  as  follows: 

continuity  ~(pu  + <pV»  + ^(pv)  0 (2.13) 

it- momentum  (pu  + <pV»~  * pv  ~ + 


- -T(p<uVW„V» 


- g--(p( u'v')tv(/iV»  »-0(€*) 


y- momentum  (pu  f (pV ))  +pv  +■ 
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3-y:(p<vV»  I0(tl)  (2.15) 

turbulent  kinetic  (/xi  r (pV»  “ rpv 
energy  y 

- |y«P V>+p<cV>  + 

(p'cVjJ-p^  + Cfc3)  (2.16) 

where  the  ( ) signify  mean  values.  Here  V Is 
written  to  Include  the  mass  flow  correlation  and  is 
defined  as 

pv*pv  + (oV)  (2.17) 

(The  bar  will  be  omitted  from  v in  all  subsequent 
equations.)  The  equation  of  state  is  written  In  the 
form 


P — PT 


The  thermal  energy  is  expressed  as  a condition  of 
constant  total  temperature  by 

t*T  <*•»> 

In  the  above  equations  u,  v.  p,  p,  T,  pan  I > are  re- 
spectively, the  poiwiimensioiul  mean  streamwise 
and  normal  component  of  velocity,  density,  pres- 
sure, temperature,  viscosity  coctficier.t,  and  ra- 
tio of  specific  heats.  Primes  denote  non-dimen- 
sional fluctuations.  The  quantity  c is  the  turbu- 
lent kinetic  energy  and  <d  is  the  turbulent  dissipa- 
tion function.  Velocities  are  made  non-dimen- 
sional with  re.-qicct  to  ul . density  with  pt , pres- 
sure with  plut=,  temperature  with u’VyM*  (fl  be- 
ing the  gas  constant),  turbulent  kinetic  energy 
with  u!1  and  dissipation  with  ulVlA  The  coordi- 


nates  k and  Jr  in  the  aliove  equations  at  non-di- 
mensionulized  with  respect  tc  L#. 

The  turbulent  kinetic* energy  equation  explicit- 
ly includes  certain  mean  dilatation  terms  in  a 
form  suggested  by  Bradshaw. tMl,<rn  The  factor 
«,  multiplying  these  terms  in  eq.  (2. 16)  is  a func- 
tion of  the  normal  stress  correlations  <u  V), 

(vV)  and  the  density-velocity  correlations  <p'u\ 
0>V)<  Bradshaw  has  shown,  from  a considera- 
tion of  the  magnitude  of  these  terms,  that  o 
should  be  a constant  tliat  is  approximately  - 10/3. 
It  follows  that  the  above  terms  will  first  effect  the 
present  solution  in  the  second-order  terms  of  the 
blending  layer.  This  will  have  a direct  bearing 
on  the  second-order  skin  friction  solution.  How- 
ever, as  noted  in  Refs.  26  and  27,  the  effect  of 
the  mean  dilatation  term  in  compressible  flows 
with  steep  pressure  gradients  can  be  numerically 
much  larger  than  their  formal  order  of  magnitude 
would  suggest.  Bradshaw  has  also  indicated  that 
these  effects  can  be  modeled  as  "extra  production 
terms"  through  an  increase  in  or.  These  consid- 
erations directly  apply  to  the  present  problem  and 
suggest  that  careful  consideration  of  the  choice  of 
a will  be  required  to  obtain  accurate  skin-friction 
solutions. 

The  form  of  the  thermal  energy  equation  given 
in  eq.  (2. 19)  is  a statement  that  the  total  temper- 
ature Is  constant.  The  contritwtion  from  the  v 
component  of  velocity  to  the  total  temperature  is 
neglected,  as  it  does  not  contribute  to  any  order 
considered  in  the  three  regions.  This  form  of  the 
energy  equal  ion  is  a valid  approximation  for  adia- 
batic walls  and  for  laniiii  h and  turbulent  Prandtl 
numbci n both  equal  to  one  Imt  Ref.  28*).  This 
appioxiniution  is  not  essential  and  the  pieaent 
theory  can  be  generalized  to  include  a lull  treat- 
ment of  the  energy  equation. 

The  above  system  of  equations  still  tin  ludes  a 
rather  large  number  ol  imlai'Wit  rot  relation  func- 
tions and  hence  does  i,ui  (mm  a closed  system  of 
equations.  However,  must  of  these  tmjcl.dimis 
will  drop  uni  of  the  1 nv*i  older  i pi.it  iors  in  the 
three  layei:  in  tin*  limit  ol  large  Reynolds  number 
(using  the  fluctuation  estimates  given  previously). 
The  only  survivin',  corrida l ions  apt>eur  in  tac  scc- 
ond-oiUcr  equations  for  the  blending  layer  and 
these  aie  the  ones  tli.irut.sed  by  Bradshaw, 
lu  this  region,  tire  pressure  diffusion  term,  as 
modeled  by  Bradsiiaw,  is  of  higher  order  and  ran 
be  neglected.  Also  inici  initlency  corvee  turn* 
need  not  lie  considered.  Tlie.se  factors  signifi- 


cantly simplify  the  closure  oblem  and  make 
strong  interaction  problems  more  tractable  than 
standard  boundary-layer  problems  (from  a clo- 
sure point  of  view).  The  specific  closure  approx- 
imation for  the  few  remaining  terms  will  be  given 
In  Section  V. 

In  the  following  sections  we  will  expand  the  so- 
lution in  each  region  In  powers  of  < and  logc.  It 
is  important  to  note  that  c is  related  to  the  Rey- 
nolds number  through  the  skin  friction  solution  in 
the  approaching  boundary  layer.  It  follows  that  < 
is  related  to  ReL  by 

c«0(l/logRet) 

Thus,  the  expansions,  in  terms  of  Reynolds  num- 
ber, proceed  in  powers  of  (l/logReL)and  log(l/ 
log  ReL).  As  a result,  the  approach  to  a large 
Reynolds  number  limit  could  be  very  slow.  This 
is  undoubtably  related  to  the  difficulty  in  estab- 
lishing a clear  Reynolds  number  trend  in  the  ex- 
perimental investigation  of  strong  interaction 
problems.  ' * 

To  establish  a frame  of  reference,  we  note 
that  for  zero  pressure  gradient  and  transonic 
Mach  numbers,  c is  in  the  range  0.02-0.04  for 
RcL  - 10*~  107.  Hence  c is  small  enough  (though 
not  impressively  small)  for  an  asymptotic  result 
to  be  useful. 

HI.  Initial  Profile 
(non-interacting  boundary  layer) 

One  would  expect  that  the  formal  asymptotic 
(booties  developed  for  incompressible  turbulent 
flowsm,-,,*’>  could  be  extended  to  the  cninprcssi- 
ble  boundary  layer  upstream  of  lire  interaction. 
This  solution  would  then  provide  initial  conditions 
Hu  t arc  asymptotically  consistent  with  the  math- 
ematical description  used  to  solve  the  intci action 
problem.  Afzal‘w*  has  considered  the  extension 
of  YajnikV16'  incompressible  solution  to  com- 
pressible turbulent  boundary  layers.  He  has 
shown  that  Hu*  structure  of  the  solution  is  th<> 
same  as  for  incompressible  flow.  Two  layers 
were  required,  an  outer  defect  region  and  an  in- 
ner wall  layer  with  the  same  length  scales  as  the 
incompressible  case.  Following  Yajnik’s  proce- 
dure fur  incompressible  flow,  the  solution  was 
repi csmited  by  limit-function  expansions  for 
small  friction  velocity  in  each  region.  It  was 
demonstrated  that  the  outer  expansion  of  the  inner 
solution  had  the  same  general  form  as  the  inner 


"Analysis  of  Uotta’fl  work  indicates  that  the  effect  of  the  interaction  docs  not  change  this  conclusion 


expansion  of  the  outer  solution.  This  indicated 
that  there  r/as  the  possibility  of  term  by  term 
matching  of  the  two  expansions.  However,  Afzal 
worked  with  an  unclosed  system  of  turbulent  equa- 
tions and  did  not  show  that  the  two  expansions  did 
indeed  match.  In  the  present  study,  we  were  able 
to  show  that  the  expansions  employed  by  Afzal  did 
not  match.  The  difficulty  stems  front  the  fact  that 
the  density  varies  by  0(1)  from  wall  to  free 
stream  values  across  the  wall  layer.  In  the  lim- 
it-function expansions  of  Rof.  19,  the  density  is 
treated  us  a constant,  equal  to  its  wall  value,  in 
the  lowest-order  inner  solution.  Our  results  in- 
dicate that,  although  this  Is  a valid  asymptotic  re- 
sult, it  is  not  valid  in  an  inner  region  large 
enough  to  overlap  the  outer  region.  * 

This  difficulty  can  bo  overcome  by  using  a 
generalized  expansion,  which  is  not  a limit-func- 
tion expansion,  to  represent  the  solution  in  the 
watt  layer.  To  obtain  an  overlap,  wc  retain  the 
density  terms  to  all  orders.  Van  Driest*2**  has 
shown  how  this  can  be  carried  out.  In  his  early 
work,  he  employed  the  Prandtl  mixing-length  the- 
ory to  model  the  Reynolds  stresses,  and  com- 
puted the  density  and  temperature  from  a Croceo 
integral  representation.  Exact  solution-*  of  the 
resulting  equations,  valid  in  the  fully  turbulent 
part  of  the  wall  layer  were  obtained  for  both  insu- 
lated and  non-ad iabatic  walls. 

The  Prandtl  mixing-length  equations  can  be 
shown  to  be  an  asymptotically  consistent  result, 
fqllowing  from  the  behavior  of  the  turbulent  kinet- 
ic-energy equation  in  the  fully  turbulent  part  of 
the  wall  layer.  Since  our  assumption  of  coastant 
total  temperature  is  contained  within  the  Crocco 
integral  representation,.  Van  Driest’s  solution  can 
be  made  the  basis  of  a formal  asymptotic  theory 
for  compressible  wall  layers.  We  follow  this 
procedure  in  developing  the  wr.il  layer  solutions  in 
the  interaction  region  (see  Section  VI).  However, 
in  establishing  the  Initial  profile,  ?e  found  it 
more  convenient  to  use  the  Mais'.-MacDonald*15* 
(also  see  Ref.  31)  correlations  of  compressible 
boundary  layer  profiles.  They  employed  the  Van 
Driest  solution  in  the  form  of  a transformation 
from  compressible  to  incompressible  profiles. 
Van  Dricst’s  generalized  solution  for  an  insulated 
plate  can  be  wr’Uen  as 


u =-sinaa 
a 


a‘-if ^ ■tLm-  «•» 

where  u is  the  compressible  velocity  and  u is  the 
related  incompressible  profile.  The  density  and 
tempeudure  can  be  obtained  from  the  total  tem- 
perature (constant  in  the  present  analysis)  as  fol- 
lows: 

1-aV  (3.3) 

P 1 * 

Van  Driest  *s  solution  for  u was  valid  only  in 
the  fully  turbulent  part  of  uu:  wait  layer,  Maine 
and  Macl-tonald*40'  made  a significant  contribution 
when  they  recognized  !l*;R  a lame  number  of  com- 
pressible profiles  can  be  correlated  by  identifying 
I*  with  Coles’  law  of  the  wake- law  of  the  wall  In- 
compressible correlation.  t“s*  Thus,  they  as- 
sumed that  the  compressible  profile  was  given  by 
cq.  (3,1)  with  Q from  (w.ing  a non-dimensional 
form  of  Alber  and  Coats"1*  notation) 


U -< 

(u. ; 


logy  c ~(1  * rosny)J 


Equivalently,  In  inner  variables 


u - u t ^ log  y * i-  B0  + - (1  - cosiry ) j (3.5) 


u, -^sin'‘a  (3.6) 

and  a Is  defined  by  Eq.  (3.2)  and  the  Van  Karman 
constant,  k = 0. 41 . The  quantity  uT  is  the  non-di- 
mensional friction  velocity,  based  on  wall  density 
and  is  related  to  our  small  parameter  < by 

(3-7) 

The  wall  density  can  be  found  by  evaluating  cq. 
(3.3)  in  the  free  stream, 

P.  - l - a*  (3. 8) 

The  Inner  coordinate  y*  Is  a stretched  wall  vari- 
able defined  by 


•This  conclusion  was  also  arrived  at  by  Adamson  and  Fco.  *2h  A general  discussion  of  this  type  of  match- 
ing problem  Is  found  in  an  article  by  Lagcrstrom  and  Castcn.  *M* 
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where  60  the  boundary  layer  thickness.  The 
constants  B0  and  i arc  the  usual  ones  appearing  in 
Cole's  formulation.  Thus,  B„-5.0  and  for  zero 
pressure  gradients,  i -0.5. 

The  skin  friction  law  for  compressible  flow 
(e.g.  Kef.  31)  follows  from  eqs.  (3.4)  and  (3,5) 

as 

'u.  1,  r/V^K.  I „ 25  „ in. 

*~I*L_5rTB,t7  (,-10) 

where  R#0  is  the  Reynolds  number  based  on 
boundary- layer  thickness  and  free  stream  quanti- 
ties. The  Maise- McDonald  representation  can  be 
viewed  as  a two-parameter  fit  of  the  upstream  ve-. 
locity  profile.  The  two  parameters  can  be  our 
•kin  friction  parameter  c and  the  Colt.  ' shape  pa- 
rameter 5.  The  skin- friction  law  then  provides  an 
expression  for  the  boundary- layer  thickness  in 
terms  of  c and  5.  Namely, 

”•» = <t£:cxp{< [“• ■ (Bn *t)]|  (3-n) 

The  Maise-McDonald  profile  ran  be  written  in 
the  form  of  a limit-function  expa  ion  in  the  de- 
fect layer  by  substitute  cq.  (3.-*)  into  eq.  (3. 1) 
atu  • xpanding  the  result  for  c-0  keeping y fixed. 
This  yields 


“o«ttr  = l ♦cu,(y)4c*~2,u5(y)  ♦ ••• 


1 JT 

u,(y)-  - lop.y--  (l  *comry)  . 


The  density  and  temperature  profiles  must  be 
computed  in  the  inner  layer  from  the  exact  ex- 
pression given  in  eq.  (3.3),  in  order  that  the  in- 
ner and  outer  expansions  overlap. 

Wc  call  attention  to  the  fact  that  the  above  ex- 
pansions are  not  valid  in  the  viscous  part  of  the 
wall  layer.  To  obtain  initial  profiles  in  this  re- 
gion, we  would  have  to  extend  Van  Driest’s  equa- 
tions to  Include  viscous  terms.  Although  this  ap- 
pears feasible,  it  would  involve  numerical  inte- 
grations which  add  considerable  complexity  to  the 
solution.  Fortunately,  this  part  of  the  initial  pro- 
file will  not  be  required  in  the  present  investiga- 
tion. 

IV.  Outer  Layer 

Details  of  the  flow  in  the  outer  or  defect  re- 
gion (i.e.  se$  2)  will  °°w  be  considered. 
Stretched  variables,  x,  y are  introduced  using  the 
length  scales  defined  by  eqs.  (2.9a,b) 

x*  (4.1) 


where  is  the  (dimensional)  thickness  of  the  ap- 
proaching boundary  layer.  It  is  the  only  length 
entering  into  the  problem  and  can  be  determined 
from  the  velocity  profile  upstream  of  the  interac- 
tion. We  introduce  these  length  scales  into  the 
governing  equations  given  by  eqs.  (2..13)  through 
(2.19)  and  assume  Jthc  non* dimensional  solution  in 
this  region  can  be  expanded  in  the  form 


1 

In  Appendix  A,  we  demonstrate  that  expansion 

u 

I c ni|(x,y  * ... 

(4.3a) 

Y 

of  the  Maisc-M*  )Xm::ld  profile  in  a hnot-hiudion 
form  in  the  wall  layci  (i.e.  c-0,  y*  fixed)  lesnlts 

V 

c^VjU.y)*  ... 

(4.3b) 

I) 

In  an  inner  solution  that  does  not  mutch  to  the  out- 
er solution  given  in  cq.  (3. 12).  Wc  also  show  that 

P-1 

I*.  ♦ <Pi(x,y) *»  ... 

(4.3c) 

the  generalized  inrn*i  solution  for  u,  lormed  by 
(In  substitution  of  eq.  (3.M  into  :q.  (3.1),  does 

P 

1 • rp,(\,y)  i . . . 

(4.3d) 

mulch  to  cq.  (2.121.  Thus,  acceptable  (t.c. 
mat*  liable*)  two  layer  .solution  for  the  non-  in**-*  - 

T 

T»  i iTl(x, y)  i . . . 

(4.3e) 

I) 

, 

acting  iKiuudai  y layer  is  given  by  the  limit  fum  - 
tioii  expan  don  eq.  (3.12)  fur  the  Older  icgionaiid 

? 

€rt^(y>  ♦ 0(c)) 

(4.30 

by  the  generalized  expansion  eqs.  (J.  11)  and  (3.  !j) 
in  the  wall  layer. 

r 

<*(‘-ii(y)  * 9(c)] 

(4.3g) 

5) 

For  future  reference,  we  give  the  outer  expan- 
sion for  the  initial  density  and  lemjieraturo  pro- 

where the  quantity  ; is  the  non-dimeuslonal  Rey- 
nolds stress,  defined  by 

files; 

T - — p (u  *V 

(4.4) 

V leh- 
* \ 


p l<<p,<y)<0(ca)  1 » c 2m,u/y)  < <Kt’>  (3.14) 

T-tIfU  - <2m#ul(y)j  ♦ 0(cO  (3.15) 


The  magnitudes  of  Uk?  fim-order  terms  of  u,  p,  p 
and  T are  determined  from  the  strength  of  the 
impinging  sliuckwave;  the  order  of  v follows  from 


mass  flow  requirements  and  the  transonic  length 
scales  given  in  cqs.  (4.1)  and  (4.2).  Effects  of 
turbulent  stresses  do  not  enter  the  expansion?  for 
u,  v,  p,  p and  T until  terms  of  0(<i/2),  hence  only 
Inviscid  equations  need  be  considered  to  this  or- 
der. This  considerably  simplifies  the  problem, 
since  large  numbers  of  unknown  correlations  need 
not  be  considered  in  the  outer  region. 

The  lowest-order  Reynolds  stress  rQ(y)  and 
turbulent  kinetic  energy  e„(y>  are  determined,  in 
principle,  from  Ihe  upstream  profile.  Fortunate- 
ly these  terms  are  not  required  to  determine  the 
leading  terms  of  the  outer  solution.  Considera- 
tion of  the  advection  terms  i. • the  turbulent  kinet- 
ic-energy equation  demonstrates  that  the  turbu- 
lent stresses  are  frozen  at  their  upstream  values 
to  0(c)  as  Indicated  iu  cqs.  (4. 3f)  and  (4. 3g).  In 
the  present  investigation,  only  the  leading  terms 
in  the  outer  Inviscid  (low  will  be  considered, 

Substitution  of  cqs.  (-1. 1)  through  (1. 3)  into 
eqs.  (2. 13) through  (2. 19)  and  carrying  out  the 
limit  c-0,  Xt  fixed  leads  to  the  following  set  of 
equations  governing  the  first-order  solution: 


(Xi  - (y  - Du,  ♦ p%  ♦ (y  ♦ l)u,Ju,t  - v,r  - 0 
du. 

%-v>.  dif 

p1=ul(y)-ul 

Pi=Pt(y)+Pi 

Tj  *-  (y-  Du,(y)  + (y-  Dp, 


(4.5) 

(4.6) 

(4.7) 
(4.0) 
(4.9) 


where  the  initial  density  distribution  p,(y)  is  given 
byeq.  (3.14)  as 


p»(y)»(y-l) ut(y), 


(4.10) 


U|(y)  is  the  defect  part  of  the  initial  profile  given 
by  eq.  (3.13)  mid  x»  is  a transonic  turbulent  inter- 
action parameter  defined  in  eq.  (2.6),  repeated 
below  as 


Mi-1  Mi-1 

*~r— 


(4.11) 


Notice  that  the  outer  solution  exhibits  a funda- 
mental similarity  in  that  it  only  depends  upon  Uu.- 
two  basic  parameters  c and  M„  In  the  combination 
appearing  in  the  definition  of  Xt  • This  indicates 
that  the  solution  is  similar  for  fixed  values  of  xt . 

Equation  (4.6)  indicates  that  a potential  func- 
tion can  be  introduced  and  the  solution  for  the  ve- 


locity components  expressed  the  form 

u,  ^u,(y)  + 0a  (4.12) 

v,=*,  (4.13) 

Thus,  the  solution  in  the  defect  region  is  an  irro- 
tationul  perturbation  to  a weakly  sheared  oncom- 
ing (rotational)  stream.  The  potential  function 
satisfies  the  following  generalized  transonic  flow 
equations: 

(yiDIMy)-^)^  I***0  (4.14) 

where 

x(y)"“["Vuij  (4-15) 

The  boundary  conditions  for  eq.  (4.14)  can  be 
found  by  matching  to  the  solution  iu  neighboring 
regions.  For  large  y,  the  solution  must  approach 
the  discontinuous  normal  shock  solution  which  in 
the  prcscut  notation  leads  to  the  condition 


4>*:i 


0 x<  0 , y-  «o 
h\m  x>0,  y ~ 


(4. 16) 


with  A.  - - x*  /y  *■  1 and  where  x -0  is  the  shock  lo- 
cation in  the  non-interacting  inviscid  flow.  Since 
the  perturbation  induced  by  the  interaction  must 
vanish  far  upstream  and  downstream  of  the  Im- 
pingement point,  we  have 


Urn  <Mx,y)*0 
lim0,(x,y)»2X. 


(4.17) 

(4.18) 


The  formulation  of  the  boundary -value  problem  Is 
completed  by  the  specification  of  surface  condi- 
tions at  y - 0.  This  is  obtained  by  matching  to  the 
Inner -layer  solutions  in  the  following  two  sec- 
tions. Anticipating  the  results  of  these  sections, 
we  have 


$,(x,  0)  = 0 


(4.19) 


Eq.  (4. 19)  implies  that  displacement  effects  In- 
duced by  the  inner  layers  do  not  contribute  to  the 
v velocity  to  0 (t,/2).  This  will  be  confirmed  in 
the  following  sections. 

The  basic  theory,  defined  by  cqs.  (4. 14)  and 
(4. 15)  is  useful  over  a limited  range  of  hlarh 
numbers  near  one  (e.g. , M.  ^ 1.2).  The  accu- 
racy of  t!»e  theory  can  be  considerably  Improved 
by  retaining  the  exact  Mach  number  dependence  of 
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a number  of  coefficients  in  the  expansion  that, 
formally,  should  be  set  equal  to  these  limiting 
values  at  M.=  l.  Specifically,  the  outer  layer 
governing  equation  is  generalized  to  the  following: 

[xl+(y+i)Mlul+(y+i)Ml«>I]<>.,-^nr,!0^  ^ 

The  boundary  conditions  for  eq.  (4.20)  are  un- 
changed except  for  the  far -field  normal  shock 
condition  given  by  eq.  (4. 16).  The  jump  condition 
associated  with  eq.  (4. 20)  is 

( 0 , x<0 , y-« 

"(irSto.  -• 

Which  contains  an  extra  Mach  number  term.  Al- 
though the  Mach  numbers  appearing  explicitly  in 
eqs.  (4.20)  and  (4.21)  could  be  formally  set  equal 
to  one,  their  retention  will  lead  to  more  accurate 
solutions.  The  original  equations  given  in  eqs. 
(4.14)  and  (4.15)  will  be  referred  to  as  the  basic 
theory  while  the  new  eq'iations  will  be  referred  to 
as  the  extended  theory.  Improved  accuracy 
arises  because  the  new  jump  conditions  yield  the 
exact  value  of  the  static  pressure  rise  behind  a 
normal  shock  wave.  An  important  feature  of  the 
extended  theory  is  that  solutions  of  the  new  equa- 
tion can  be  obtained  from  a simple  transformation 
of  solutions  to  the  basic  small -disturbance  equa- 
tion. In  Sec.  VII,’  we  develop  a numerical  tech- 
nique for  solving  the  basic  small -disturbance 
equations  given,  iu  eqs.  (4.  14)  and  (4.15).  These 
cart  be  converted  into  solutions  of  the  more  accu- 
rate extended  equations  through  the  following 
transformation.  If 

p,  = Pi(x,y;Kt)  (4’2Z) 

is  a solution  for  the  pressure  obtained  from  top- 
(4.14)  and  (4.15)  for  a given  xt  equal  toK»,  then 
the  solution  of  the  extended  equations  is  given  by 

pr.^Pi(M-x,y;K,)  t4’23^ 

where  Kt  is  related  to  the  Mach  number  by 


Kt=  "cmT 


Cf(x*,y*;c,M„,  t,r) = cCt(s,  y;K»,  *)(4. 25) 
where  s and  Kt  are  defined  by  the  relations 

£.ct',*«c,/iii „(y+l)s  (4.26) 

6J 

_ Mi  -1  U Ml 

ani  where  C,  is*  a universal  function  of  s and  y for 
given  vaHies  ofKt  and  ».  It  follows  that  s and  y 
are  the  basic  non-dimensional  coordinates  for  the 
problem  andK,  Is  a general  similarity  variable 
for  transonic  turbulent  interactions.  The  above 
similarity  applies  for  the  law  of  the  wall-law  of 
the  wake  form  assumed  for  the  Initial  velocity 
profile.  The  similarity  holds  for  general  initial 
profiles  provided  the  flow  has  the  same  Initial  pro- 
file as  a function  of  y*/6J. 

The  present  formulation  ic  valid  for  shock  _ 
waves  of  moderate  strength,  that  is,  for  X* -«(*)• 

It  is  of  interest  to  inquire  into  the  behavior  of  the 
solutions  in  the  limits  of  large  and  small  values  of 
X . For  x««0(l),  the  sonic  line  is  in  the  main 
part  of  the  defect  layer  and  the  appropriate  length 
scales  are  the  ones  used  in  the  present  analysis. 
The  sonic  line  approaches  the  wall  for  increasing 
values  of  x,  and  moves  into  the  inner  layers  for 
large  xt*  In  the  large  xt  limit,  it  can  be  shown 
-that  the 'flow  In  the  defect  layer  is  described  by  a 
■ linearization  ‘about  the  piecewise  constant  normal- 
shock  solution.  A boundary-value  problem  simi- 
lar to  that  treated  by  Adams"”  arises.  This  solu- 
tion must  develop  a singularity  at  the  wall.  Ihis 
suggests  that  the  solution  in  the  inner  layer  for 
large  x,  must  be  relatively  complex  and  probably 
involves  normal  pressure  gradients-  An  analysis 
of  this  problem  has  not  yet  been  carried  out. 

For  the  weak  shock  limit  of  xt-0»  the  present 
formulation  should  reduce  to  that  of  Adamson  and 
Feo.  ttl)  This  limit  process,  however,  involves  a 
complicated  singular  perturbation  problem  which 
) is  considered  in  Appendix  B.  There,  we  show 
that  the  present  analysis  includes  the  weak  shock 
theory  of  Ref.  21. 


V.  Blending  Layer 


Then,  transformations  provide  a simple  and  dircet 
method  for  greatly  improving  the  accuracy  o the 
small-disturbance  theory.  Equation  (4. 20)  also 
leads  to  a turbulent  similarity  theory,  for • «- 
ample,  it  follows  from  eqs.  (4.20)  and  (4.21)  tt  at 
the  pressure  coefficient  cun  be  written  m the  form 


The  structure  of  the  solution  in  the  two  inner 
layers  is  such  that  the  first  two  terms  of  the  solu 
lion  in  the  wall  layer  can  be  conipleteiy  deter- 
mined without  consideration  of  the  blending -l.»> 
solution.  The  first  two  terms  in  the  skin  friction 
can  be  obtained  by  comparing  the  wall-layer  solu- 
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toon  directly  with  the  solution  in  the  outer  or  de- 
fect region.  In  this  way,  it  can  be  shown  that  the 
change  in  skin  friction  due  to  interaction  with  the 
shock  wave  is  the  order  of  the  pressure  rise, 
which,  in  the  present  problem  is  0(c).  Thus  the 
shock  strengths  considered  in  the  present  inves- 
tigation, xt  =0(1),  are  not  large  enough  to  cause 
separation. 

Because  of  the  freezing  of  the  Reynolds 
stress  in  the  outer  layer,  as  discussed  previ- 
ously, the  outer  and  wall  layer  regions  dr  not 
overlap.  Hence,  to  be  logically  consistent  and  to 
confirm  the  skin-friction  solution,  wc  should  com- 
plete the  wall-layer  solution  by  matching  to  the 
solution  in  the  blending  layer.  Therefore,  in  the 
present  section,  we  will  develop  the  asymptotic 
solution  in  the  blending  layer,  carrying  all  terms 
that  contribute  to  the  second-order  skin  friction. 
For  greater  generality,  we  will  retain  all  Mach 
number  terms  in  these  expansions  and  also  in  the 
wall  layer.  Transonic  approximations  consistent 
with  the  outer  solution  can  then  be  obtained  by  ex- 
panding these  solutions  in  powers  of  M*. 

It  is  only  at  this  stage  of  the  analysis  that  em- 
pirical assumptions  need  be  made  in  order  to 
close  the  system  of  equations.  Most  of  the  cor- 
relations that  appear  in  the  time -averaged  com- 
pressible flow  equations  vanish  to  the  order  re- 
quired in  the  present  study.  Four  assumptions 
must  be  made  in  order  to  close  the  blending -layer 
equations  to  second  order.  These  are  U)  an  as  - 
sumption  relating  turbulent  kinetic  energy  to  the 
Reynolds  stress,  (2)  a model  for  the  dissipation 
function,  (3)  an  estimate  of  the  pressure  diffusion 
terms  and  (4)  an  estimate  of  the  mean  dilatation 
parameter  a.  In  order  to  conveniently  arrive  at 
a definite  set  of  equations,  we  will  follow  Brad- 
shaw,s<,<,2*)  model  of  these  terms.  Thus,  the 
pressure  diffusion  terms  are  asymptotically  of 
higher  order  and  do  not  appear  in  the  second -or- 
der equations.  This  seems  to  be  a reasonable  ap- 
proximation, since  experiment  al  data  indicate  that 
the  pressure  diffusion  terms  a.*e  much  smaller 
than  the  other  terms  appearing  in  the  turbulent  ki- 
netic-energy equation.  In  Bradshaw's  model,  the 
turbulent  kinetic  energy  is  assumed  to  be  given  by 


L«(y) 


($.3) 


@«*T/2a  tp 

where  his  suggested  value  of 
at*0.15 


(5.1) 


(5.2) 


The  dissipation  length,  L Jy),  is  assumed  to  have 
the  following  expansion 

<3  .4) 

where  k is  the  von  Karman  constant  (k-0.41)  and 
li,lt  are  constants  that  can  be  determined  from  a 
curve  fit  of  Bradshaw's  suggested  function.  The 
values  of  lj,i2  do  not  appear  in  the  second-order 
equations  ami  hence  are  not  actually  needed. 

The  requirement  that  the  inertia  and  Reynolds 
stress  terms  in  the  streamwise  momentum  equa- 
tions arc  of  the  same  order  leads  to  the  lengths 
given  in  eqs.  (2.0a)  an  ! (2.12).  Using  these 
length  scales  the  following  stretched  variables  are 
defined  as 


x*  = cl/%*x 


(5.5) 

(5.6) 


Preliminary  analysis  of  the  x momentum  equation 
indicates  that  the  first-order  Reynolds  stresses 
are  balanced  by  the  second-order  inertia  terms. 
This  complicates  the  development  of  the  solution, 
in  that  the  expansion  for  u and  p must  be  carried 
out  to  second  order  for  a computation  of  the  first 
term  in  the  solution  for  the  Reynolds  stresses. 

The  expansion  for  u contains  a number  of  log< 
terms  that  must  be  considered.  These  terms 
arise  from  contributions  of  the  outer  solution  and 
from  the  initial  profile.  They  can  easily  be  ex- 
tracted from  these  expansions  when  written  in 
blending  layer  variables.  Taking  these  terms  into 
account  and  after  some  simple  manipulations,  the 
expansions  for  the  blending  layer  sdution  can  be 
written: 


u = 


is  most  suitable  for  regions  near  the  wall.  The 
dissipauon  function  is  taken  in  tlie  standard  form, 


l+(^)<togc-*W)*(|J  ^ 

- ^j[m,H(Y)  - (1  + 2m,)P,(x)] 

+ <tU,(x,Y)+... 

v=t*V,(x,Y) 

P=l  + €P,(x)  + €*P,(x)  + ... 


P = 1 + y J-  €log«  +<[2m,n(Y)+MiPj(x)] 


(5.7) 

(5.8) 
(5.8) 

(5.10) 
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7) 


B) 

B) 


0) 


Tsl  (y-l)€log€+€(y-lHpi(x)’’H{Y)]  + »*» 

V2*'  (5.11) 


what  simplified  by  tlie  substitution 

U|  ^ U(x,Y)  -P,(x)  4 (1 4 2m.)H(Y)Pi(x)  + i m.H*(Y) 

(5.19) 


t = c*[l  4«rl(xl  y ) 4 0(c*log<)]  (5.12) 

e= c*[2a,  4 cEi(x,  Y)  40(c*logc)]  (5.13) 

where 

H(Y)-^  (logY  ~ 2H]  . (5.14) 

Consideration  oi  the  normal  momentum  equa- 
tion indicates  that  the  pressure  is  constant  across 
the  blending  layer  to  second  order.  Thus,  the 
functions  Px(x)  and  P,(x)  appearing  in  eq.  (5. 9) 
can  be  determined  from  the  wall  values  of  the  out- 
er solution.  Assuming  that  wall-layer  displace- 
ment effects  are  negligible  to  second  order,  (as 
will  be  verified  in  the  next  section),  we  find  from 
continuity  considerations 

In  the  transonic  limit  V,  actually  becomes  first 
i>rdcr  in  t 

Thus,  the  blending  layer  is  a parallel  stream  to 
third  ord'-r  at  transonic  speeds. 

Substitution  of  the  previous  expansion  into  the 
x momentum  and  Imbuh  nt  binotic- cneiay  equa- 
tions results  in  the  following  equations; 

t&. iv) 


r,«r(x,Y) 

(5.20) 

whereby 

u,-rT=o 

(5.21) 

rg-2a1Ur  = 

--SMr-oPito]-2ai“M- 

xY 

dPj 

dx 

(5.22) 

and  where 

* 

<r*  Mi  4 2(1 4 2m.)  . 

(5.23) 

Thus 

Lim  U(x,Y)  =0 

(5.24) 

LimlXx,  Y)  = 0 . 

(5.25) 

Homogeneous  initial  conditions  are  obtained  be- 
cause < has  been  defined  to  be  equal  to  the  exact 
value  of  the  skin  friction  at  the  initial  station. 
Boundary  conditions  arc  determined  by  matching 
to  the  defect  layer  for  Y - 00  and  to  the  wall  layer 
for  Y-  0.  Matching  to  the  defect  layer  leads  to 
the  condition 

Lim  1)(>,Y)-  0 (5.26) 

r— » 

In  the  following  section,  we  show  that  the  firs* 
order  skin  friction  can  be  determined  from  the 
outer  and  waP  -layer  solutions  independent  of  the 
blending-iayer  solution.  This  leads  to  the  follow- 
ing inner  condition: 


1 (5.1  fa) 


where  the  assumptions  given  in  cqs.  (5.1)  through 
(5  4)  were  used  to  close  the  system  of  equations. 
We  call  attention  to  the  fact  that  advection  of  tur- 
bulent kinetic  energy  is  retained  as  a leading  term 
and  that  convection  is  rel, lined  in  the  momentum 
equation.  Thus,  the  shear  stress  gradient  r is 
not  constant  in  this  region.  Hence,  the  blending 
layer  is  not  a generalized  equilibrium  layer  in  the 
sense  of  Townsends  theory00  of  turbulent  separa- 
tion. 

The  equations  (6.17)  and  (5.18)  can  be  some- 


(s.«) 

whore  1’j  is  determined  from  solution  to  the  outer 
problem  and  is  a constant  determined  from  the 
wall  layer  solution. 

The  solution  to  the  Ixnindary  value  problem  for 
U(x,y)  can  bo  shown  to  have  the  following  behavior 

for  y — 0 

U = 1 It,.(x)  - cr P,(x)][log Y - 2ff]  4 C2(x) 

(5.28) 

where  C4x)  is  a function  that  can  be  extracted 
from  the  numerical  solution  of  the  above  problem. 
The  function  C2(x)  is  required  in  the  wall -layer 
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solution  to  compute  the  second -order  skin  fric- 
tion. Thus,  the  flow  in  the  blending  layer  first 
affects  the  solution  for  the  skin  friction  to  second 
order.  It  follows  that  the  effects  of  empirical 
closure  approximation  also  enter  only  'it  this  or- 
der. 

Before  leaving  this  section,  we  note  that  a 
compact  statement  of  the  blending-layer  problem 
can  be  obtained  by  Introducing  a generalized  po- 
tential, defined  by 

U=0T  (5.29) 

r*fl,  . (5.30) 

Equation  (5. 21)  is  automatically  satisfied,  and 
substitution  into  eq.  (5. 22)  leads  to  the  following 
second-order  partial  differential  equation: 


wt^re  D is  the  damping  factor.  There  have  been 
t nun. be r of  choices  for  D,  suggested  in  the  lit- 
e.-atu  .•©.  They  are  ail  equivalent  If  pressure  gra- 
dients are  not  important  in  the  wall  layer,  be- 
cause then  they  arc  all  based  on  the  same  set  of 
zero  pressure -gradient  data.  Although  the  pres- 
sure gradients  in  the  present  interaction  problem 
are  large,  the  shear  stress  gradients  near  the 
wall  are  even  larger.  An  order  of  magnitude 
analysis  of  the  momentum  equation  using  the 
present  length  scales  shows  that  the  pressure 
gradient  and  inertia  terms  are  exponentially  small 
compared  to  the  shear  stress  terms  and  hence  can 
be  neglected  to  all  orders*.  It  follows  that  any  of 
the  standard  zero  pressure-gradient  damping  fac- 
tors should  be  adequate  for  the  present  work.  For 
example,  Van  Driest  suggested  the  form 

D = 1 -exp(-yY26. ) (0.1) 


n„  - 2a, n„  n,=  ?l&(-  - 2.,,  a ^ m *.  . 


kY 


dx 


(5.31) 


Equation  (5.31)  is  a hyperbolic  equation  as 
must  arise  when  the  pressure  diffusion  terms  are 
neglected  (or  approximated  by  Bradshaw’s  model). 
The  last  term  in  eqs.  (5. 18)  and  (5.31)  arise  from 
mean  dilatation  terms  in  the  turbulent  kinetic -en- 
ergy equation.  As  noted  by  Bradshaw and 
discussed  previously,  the  values  of  a appearing 
in  these  terms  must  be  carefully  chosen  to  model 
the  effects  of  extra  mean  dilatation.  These  terms 
are  believed  to  have  a value  of  about  ten  (Ref.  26) 
and  therefore  can  be  expected  to  have  a larger  ef- 
fect mi  the  skin  friction  than  their  formal  magni- 
tude (i.e.,  c*)  would  suggest.  Considerations  for 
determining  values  of  a are  given  in  Refs.  26  and 
27. 


VL  Wall  Layer 

In  the  wall  layer,  all  correlations  except  the 
dissipation  function  are  asymptotically  negligible 
for  the  insulated  flat  plate  under  consideration. 

In  our  solution,  we  model  the  dissipation  function 
in  the  usual  manner  by  introducing  a damping  fac- 
tor to  account  for  viscous  effects  near  the  wall. 
Thus,  we  assume  the  dissipation  is  written  in  the 
form 

4 KyD 


whore  y*  is  the  local  wall  coordinate  defined  by 


y* 


J-jiaaJL 

lP:(xK(x)Jy 


(6.2) 


The  only  property  of  D actually  required  in  the 
present  work  is 


Um  D=  1 + exponentially  small  terms  (8. 3) 


The  solution  in  the  wall  layer  will  be  carried 
out  in  stretched  coordinates  x,  y*.  We  also  In- 
troduce a wall  layer  velocity  u*  by 

u»uT(x)u*  (6.4) 

We  recall  the  definition  of  the  local  friction  veloc- 
ity 


where  p„(x,c)  and  tw(x»«)  are  the  local  wall  den- 
sity and  local  skin  friction,  respectively.  Substi- 
tution of  the  above  wall -layer  variables  into  the 
governing  equations  and  carrying  out  the  limit  < 

- 0 (however,  retaining  all  Mach  number  terms) 
and  noting  the  order  of  magnitude  estimate  in  eq. 
(6.5)  we  obtain 


T 

% 


1+0(0 


(6.6) 


♦As  is  well  known06’* w6>  this  conclusion  does  not  hold  near  separation.  Hence  the  present  developments 
are  likely  to  be  not  uniformly  valid  near  such  points.  An  additional  expansion,  including  pressure  gradi- 
ents, and  a small  length  scale  is  probably  required  for  this  region. 
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wf (6-7’ 

fi«-=l-ay(u*)1  (6.8) 

P 

(6-9) 

V * O(exponentially  small)  (6. 10) 


where  K {p*/p}  Is  defined  from  a specification  of  a 
particular  viscosity  model. 

Equation  (6.6)  follows  from  the  streamwise 
momentum  equation.  It  is  the  usual  condition  that 
the  total  stress,  laminar  plus  turbulent,  is  con- 
stant across  the  wall  layer.  Equation  (6. 7)  is  a 
mixing-length  equation  for  the  Reynolds  stress, 
augumented  by  a damping  factor  to  account  for 
viscous  effects  on  the  length  scale  in  the  wall  lay- 
er. In  the  present  context  it  follows  as  a formal 
asymptotic  consequence  of  the  dominance  of  the 
production 'and  dissipation  terms  in  the  turbulent 
kinetic -energy  equation  as  y~  0.  We  note  that  the 
mean  dilatation  terms  are  negligible  to  all  orders 
in  the  wall  layer.  The  density  is  computed  from 
eq.  (6.8)  which  follows  from  the  equation  of  .state 
and  the  condition  that  the  pressure  mid  total  tem- 
perature are  constant  across  the  wall  layer.  The 
viscosity  is  computed  from  the  density  using  an 
equation  of  the  form  given  in  cq.  (6.  9). 

These  are  the  usual  equations  used  for  incom- 
pressible arid  compressible  wall  layers.  Inser- 
tion of  the  normal  momentum  equation  indicates 
(hat  pressure  is  constant  across  the  wall  layer  to 
all  orders.  Thus,  the  pressure  is  imposed  from 
the  outer  solution  ami  is  given  the  form 

P » "r]jr  ♦ «P*,(x)  ♦ c*p.f  lx)  * . . . (G.  1 1 ) 

Equation  (6. 10)  indie. lies  tin-  v eouqmnent  of 
velocity  is  <>.|W»n**ntially  small.  This  implies  that 
wall -layer  displacement  effects  are  negligible  to 
all  orders  and  that  the  noma)  velocity  linundaty 
conditions  used  in  the  outer  regions  are  correct. 

In  order  to  develop  wall -layer  expansions  ttiat 
are  valid  in  a region  large  enough  to  overlap  the 
blending  layer,  we  follow  the  procedures  adopted 
in  Section  111  for  non-interacting  boundary  layers. 
Thus,  we  transform  eqs.  (G.G)  through  (6.  8)  using 
Van  Driest’*  generalized  ve!oi  ities 

u = ufti*~  ~ siu(auru)'  (6,12) 

4 


We  also  define  new  dependent  variables  F, 
place  of  r and  p by  the  relations 

Rin 

t*tw(x;  c)F(x,  y*;  c)  = 0(c*) 

(6.13) 

p=pjx;  c)R(x,  y*;  c)=0(l) 

(6.14) 

where  rw(x;  c)  is  the  order  of  the  shear  stress  of 
the  approaching  boundary  layer  and  is  0(c*). 
Under  these  change  of  variables,  the  wall-layer 
equations  become 

F.(K(R)/R,/*]|p«l»0«’) 

(6.15) 

VF*  tKy-)«y*-^T  ♦«(<•> 

(6.16) 

1 

cos^aiyi) 

(6.17) 

with  K(R)  determined  from  a viscosity  law. 

Equations  (6. 15)  to  (6. 17)  are  to  be  solved  sub- 
ject to  the  boundary  condition 

u = 0 for  y*  = 0 (6.18) 

We  note  that  the  density  appears  only  in  the 
second  term  of  eq.  (6.15).  This  term  is  multi- 
plied by  a viscous  term  which  becomes  exponen- 
tially small  in  the  fully  turbulent  part  of  the  wall 
layer.  As  a result,  the  troublesome  density 
terms  cause  no  mathematical  difficulties  and  we 
can  represent  the  solution  of  the  transformed 
equations  in  a limit -function  expansion. 

To  solve  Urn  resulting  equations,  we  must 
choose  a particular  model  for  the  function  D(y*) 
and  for  the  viscosity  law.  The  resulting  equations 
for  a given  choice  of  D(y’)  usually  can  only  Ik*  in- 
legiated  numerically.  However,  in  order  to 
match  to  the  bli  nding  layer  and  to  determine  the 
skin  friction,  we  only  require  a solution  in  the 
tally  turbuP  ii4.  part  of  the  wall  layei.  This  is  ob- 
tain* d in  the  usual  fashion  by  dropping  the  viscous 
term  in  eq.  (6.1b),  .setting  ’)-l.  The  resulting 
equations  reduce  to  the  ianiiliar  incompressible 
form  which  can  be  integrated  in  rlosed  form  to 
yield 

u^logy%I3(x;  c)  (G.  19) 

where  B(x,*  c)  is  a function  that  can  only  be  deter- 
mined from  an  integration  of  the  full  equations 
(i.e.,  cqs.  (6. 15)  through  (G.  17)).  Thus,  P(x;  c) 
must  depend  upon  the  density  variations  through 
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the  second  terra  in  eq.  (6, 15).  Noting  that  u,  * 0(c), 
eq.  (6. 1?)  indicates  that  D must  be  of  the  form 

B*Of+ctB,(x)****  (6.20) 

where  the  first  terra  can  be  identified  with  the  in- 
compressible value  and  the  second  term  contains 
the  first  effects  of  compressibility  in  the  wall  lay- 
er (in  transformed  variables).  The  incompressi- 
ble constant  B#  is  a function  of  the  damping  model 
employed  in  the  computations.  These  models  are 
usually  arranged  to  give  the  experimentally  ac- 
cepted value  of  B#=  5. 0. 

Thus,  the  solution  in  the  fully  turbulent  pari  of 
the  wall  layer  can  be  written  as 

u=isin  f V^io[« tog  y%!V  ''D*  * '•]} 


The  wall  density  can  be  related  to  the  surface 
pressure  through  the  equation  of  state  evaluated  at 
the  constant  wall  temper  . The  skin  friction 
is  an  arbitrary  function  oi  * and  c at  this  stage  of 
the  analysis.  It  will  eventually  be  determined  by 
matching  the  wall  and  blciMling -layer  solutions. 


Matching  considerations  suggest  that  tw 
can  be  expanded  as 

and  pv 

T.  3 <*U  ♦ < rml  * <2log«  r„n  f <*rW2l  + • • 

•1  (6.22) 

Pw-Awo  + ep^  + ^P.,**** 

(6.23) 

From  the  equation  of  state  we  find 

Pw0“l-a* 

(6.24) 

P»,  =vM*4>Wl(x) 

(6.25) 

P*,=yMipWj(x) 

(6. 26) 

sion  is  further  expanded  into  limit -function  form, 
it  can  be  shown  that  the  resulting  solution  will  not 
match  with  the  blending -layer  solution.  The  dem* 
onstration  that  the  generalized  wail- layer  solu- 
tion given  above  can  be  matched  term  by  term  to 
the  blending -layer  solution  is  somewhat  involved, 
and  the  details  are  given  in  Appendix  C.  Here,  we 
present  the  solution  for  the  skin  friction  to  second 
order  which  comes  from  the  requirement  that  wall 
and  blending-layer  solutions  match  to  second 
order 


-*--1  >< 
C*o 


+‘[yM*  *q]IV>  * 

<e.2») 

q = (6.20) 


The  evaluation  of  the  second -order  skin  friction 
from  the  above  expression  requires  a numerical 
solution  of  the  boundary -value  problem  formulated 
in  the  blending  layer  in  order  to  determine  the 
function  C2(x),  appearing  in  eq.  (6.28).  The  tran- 
sonic limit  of  eq.  (6. 28)  is  obtained  by  setting  M. 
to  one  in  the  above  expression.  For  reference, 
we  note 


a*0. 40825 
o~0. 94034 


| for  II.- 1. 


The  first  term  above,  is  the  wall  density  in  the 
non-interacting  flow  upstream  of  the  shock  wave. 

With  these  results  the  wall -layer  solution  for  the 
velocity  can  be  written  as 

U * a Si°  (v£o {[«  !<>«  y*  + B0]  ♦ i €ItWi  - pWJ ) ~ log  y * 

+ a<llog<(rW||)^  logy* 

* i«*( T*u  - P.„  ♦ 2p-f  - i(Tw.  - P.,)2  + 2D*)  + * • * \) 

(6.27) 

Notice  that  the  solution  for  u is  not  in  the  form  of 
a limit -function  expansion.  If  Uic  above  expres- 


VH.  Num  irical  Analysis  of  the  Cuter-Layer 
Equations 


This  section  presents  the  numerical  analysis  of 
the  first-order  (inc)  outer-layer  equations  and 
boundary  conditions  as  described  by  eqs.  (4. 14)- 
(4. 19).  As  shown  in  Sections  V and  VI,  normal 
pressure  gradients  in  the  wail  and  blending  layers 
are  aero  to  second  order.  Thus,  the  solution  of 
the  above  outer -layer  equations  will  trace  the  de- 
velopment of  the  pressure  distribution  from  a dis- 
continuous jump  at  the  normal  shock  to  a smooth 
rise  at  the  surface.  This  surface  pressure  dis-  ■ 
tribution  can  then  be  used  to  solve  the  blending 
and  wall  layers  for  the  skin  friction. 

The  complete  boundary -value  problem  is  re- 
stated here  as 

(y  + l)(*(y) (7.1) 


«»,(»)] 


Inserting  *&.  (7.  7)  and  (7. 8)  into  (7. 1),  retain* 
ing  lowest -order  terms  and  comparing  similar 
terms  in  y gives 


C|(x)  = - bo  0«*k» 


, . ( 0 , x < C 

Um«,(x,y)  = |2^x>( 


Lim  0,(x,y)  = O 


c#(x)  = - [b|  - |b0  - 0,(x,  0)1  $u(x,  0) . 


Lim  0t(x,y)  = 2A,, 


and  at  the  surface 


0,(x,  0)*0 


Now,  eq.  (7.8)  may  be  solved  for  0(x,O)  in  terms 
of  the  solution  at  a small  distance  y away  from  the 
wall,  with  the  terms  ct  and  c,  obtained  from  eq. 

(7. 9)  involving  x-derivatives  of  the  wall  solution 
from  a previous  iterate. 


x.*-Xt/v+1 


The  only  quantities  remaining  to  complete  the 
boundary-value  formulation  is  the  oncoming 
boundary-layer  defect  profile  u4(y).  In  principle, 
this  profile  may  be  obtained  from  experimental 
data.  However,  consistent  with  the  formalism 
used  to  develop  the  analysis  in  the  inner  layers, 
we  describe  u,(y)  as  the  non-interacting  boundary 
layer  profile  given  by  eq.  (3. 13).  Hence,  the 
boundary -value  problem  as  described  by  eqs. 

(7. 1 )— (7.  C)  will  possess  a two-parameter  family 
of  solutions,  depending  on  values  of  the  interac- 
tion parameter  x»  and  the  initial  profile  shape  fac- 
tor * . 


The  above  boundary-value  representation  is 
still  not  in  a form  suitable  for  numerical  integra- 
tion. The  tangential  boundary  conditions  at  the 
shock,  eq.  (7-  3)  must  be  replaced  by  the  appro-^ 
priate  Dirichlet  condition.  Without  loss  in  gen- 
erality, this  is  accomplished  by 


. ( 0 , x<0 

IWtx.y).  J2>_X|  x >0 


But.  eq.  (7.10)  introduces  another  numerical  dif- 
ficulty. The  limit  of  0(x,y)  = 2x.x  as  y-«>  cannot 
be  applied  as  x-  *.  We  have  chosen  to  solve  this 
problem  by  separating  the  infinite  domain  into  two 
parts.  Part  1,  x<0,  we  have 


Llm<>|(x,y)~0 


One  difficulty  will,  the  a»n#ve  formulation  lies 
In  the  implementation  of  the  surface  boundary 
condition,  eq.  (7.0.  Standard  numerical  tech- 
nique* evaluate  Ibis  condition  through  a reflec- 
tion plane  ptocidure,  *huh  requires  the  govern- 
ing eq.  (7. 1)  to  W evaluated  at  the  wall  y 0. 
However,  due  to  tin.  logarithmic  behavior  of  the 
initial  prolib  u,(y),  t Ik- coefficient  x(y>  becomes 
singular  for  y-  0. 


Lim*»(x,y)-0,  0j,(x(O)  0 


And,  for  part  2,  x>0,  we  have 


0ci.$~2A„x,  LinW'*(x,y)-- 0 


Limtf*(x,y)  0, 


One  way  of  circumventing  this  problem  Is  to 
expand  x(y)  foi  small  y as 


The  interface  equations  across  the  boundary  x - 0 
are 


A(y)  b0  logy  i bj  f0(y*) 


$i(0,y)"<MO,y) 


where  bQ  and  b,  are  known  constants.  The  behav- 
ior of  <>  near  the  wall  can  be  seen  to  be  of  the 
form 


<*L(0fy)- £it(0,  y)*2X. 


0(x, y)  • 0(x,  0)  ♦ c.jOOy2 logy  4 c^xly*  4 0(y*) . 


The  equation  governing  % is  identical  to  equa- 
tion (7. 1)  replacing  <>  by  <«, . The  equation  for  <*, 


(7. 8)  is  written  as 


\ 

j 


(y  * l)[My)  - 2X„  - 0J0*.  *■  $tTt ~ 0 <7. 14) 
This  formulation  is  illustrated  in  Fig.  3. 


tion  of  a discontinuous  wave  near  a caustic.  0i* 

This  latter  problem  results  in  a somewhat  similar 
equation  to  cq.  (7. 1)  with  X being  proportional  to  y. 


In  order  to  develop  a domain  suitable  for  the 
numerical  solution  of  the  above  equations,  we 
analytically  transform  the  upper  half  plane  (phys- 
ical domain)  to  a unit  square  (computational  do- 
The  particular  functions  are  chosen  to  be 

^tAjlog^log^  *A,4a*A4,  i<ik 

M*A,,  «**«**, 

^Aj  *■  Aj  log  * A4^b  ♦ A„  t>t,a 

(7.15)  * 


l A,  n,  0‘  >,v»U 

flA,  * A,  log  M!  jzfA  + At  ’Ja  » * >l Ha 

(7. 16) 

above  functions  are  seen  to  have  the  ap- 
propriate behavior  since  4(-  <*>)  = 0,  4(+«)  = 0, 
if  (O’  = 0 and  »](*)  = 1.  Furthermore,  constant  rec- 
tangular mesh  spacings  A4,  Aq  in  the  computa- 
tional domaiE:  yiela  the  variable  mesh  in  the  phys- 
ical plane  as  depleted  in  Fig.  4.  The  stretching 
functions  (7. 15)  and  (7. 16)  are  chosen  to  have  a 
uniform  grid  in  the  physical  plane  within  a rec- 
tangular region  surrounding  the  shock,  near  the 
surface.  The  mesh  lines  then  expand  nonuniform - 
iy  to  infinity  in  all  directions.  The  parameters 
Aj  - A»,  <A,  tB  and  1a  are  used  to  determine  the 
size  of  the  mesh  in  the  inner  rectangular  region 
and  to  control  the  rate  of  this  mesh  /rowth  to  in- 
finity. The  transformations  were  chosen  with  this 
uniform  inner  rectangular  grid  in  order  to  care- 
fully study  the  region  where  the  shock  meets  the 
outer  edge  of  the  boundary  layer  with  minimal  ef- 
fects of  variable  mesh  spacing.  The  derivatives 
of  the  transformation  are  used  in  the  solution  of 
cqs.  (7.1)  and  (7. 14)  as 


main). 


and 


The 


A line-relaxation  algorithm  is  developed  to 
solve  the  outer  problem.  A uniform  rectangular 
grid  is  employed  in  the  computational  plane  with 
mesh  spacings  Aij  ami  A$.  Subscripts  N and  M 
in  the  following  analysis  refer  to  the  counters  re- 
spectively locating  the  4 and  rj  mesh  lines.  For 
regions  of  subsonic  flow,  the  following  difference 
expulsions  are  used 

«■- <7-  to) 


0<r- 


se 

(7.19) 


where  the  superscript  v refers  to  values  at  the 
previous  iterate  and  v + 1 refers  to  the  present 
iterate.  The  quantity  u>  is  the  relaxation  factor 
which  on  the  basis  of  linearized  stability  consid- 
erations lies  between  0 and  2.  For  this  study, 
we  take  >1  corresponding  to  over-relaxation. 


In  regions  of  supersonic  flow,  we  take 


01s 


(7.20) 


and  for  both  regions  we  take 


6 , £LkiLZ  0Jj,Url 

2 All 


(7.21) 


(7.22) 

(7.23) 


In  order  to  determine  whether  the  flow  Is 
supersonic,  the  following  test  function  Is  used 


(7.17) 

0T?s0|^**  + 0«^«  0rr30«tlrr + 0w,4 

It  has  recently  been  found  that  an  efficient  and 
accurate  means  of  solving  this  system  consists  of 
a mixed  flow  line-relaxation  procedure  introduced 
by  Murman  and  Cole.ao>  This  technique  has  been 
successfully  applied  to  a variety  of  inviscid  tran- 
sonic aerodynamic  problems  and  to  the  cotTiputa- 


ir  i UlMf'Ti’  iriOaMtflMHl  if 


«.24) 

and  the  supersonic  differencing  procedure  is 
used  when 

Mil..- My)  >0  (7.25) 

This  particular  test  function  has  been  found  to  be 
stable  for  purely  decelerating  flows.  In  addition, 
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we  have  found  that  for  the  present  problem,  this 
choice  of  test  function  results  in  shocks  smeared 
over  a minimum  number  of  mesh  intervals  (i.e. 
1*2). 

The  above  difference  formulae,  (7. 18)— (7. 24), 
are  arranged  using  standard  methods  in  a con* 
venient  algorithm  to  update  values  of  $ in  sue  - . 
cessive  columns  starting  with  £ = 0 proceeding  to 
4*1.  Note  that  with  the  above  schemes,  the  dif- 
ference equations  will  be  non-linear  in  supersonic 
regions  due  to  the  presence  of  terms. 

This  requires  a local  fixed -point  iteration  as  dis- 
cussed in  Ref.  20. 

It  is  seen  from  eqs.  (7. 4)  and  (7. 5}  that  the 
value  of  $t-2X.*-2xt/y  + l as  x-  +«.  In  order 
to  compare  runs  at  different  values  of  Xi»  we 
found  it  convenient  to  introduce  following  addition- 
al transformations: 


(7.26) 

and 

(7.27) 

whereby  equations  (7. 1)  and  (7.2)  become 

i»(y>-$ltfir-4„  = 0 (7.28) 


and 


X(y)  = -^  + -^^u,(y)j. 

Now,  we  see  that  - - 2 as  x-  •>.  Note  that  in 
terms  of  physical  variables 

= <l/lx  = / cjf,  x v'Mi ~ l x (7 . 29) 


We  note  for  reference,  that  oq.  (7. 26)  pos- 
sesses rr  ii  eh:»racteri«lies  upstream  of  the  shock 
wave  and  altove  the  sonic  line  ot  the  em  oinim*. 
boundary  layer.  The  charaeti  ris’it  direetioiis  are 


?r-r 


(7. 30) 


We  also  define  a normalized  perturbation  pressure 
P“-  = (7.31) 

\i  X| 

which  will  be  used  in  the  next  section  where  the 
numerical  results  will  lie  discussed. 


VIII.  Results 


The  results  of  the  numerical  computation  of  ' 
the  first-order,  outer-layer  equation  are  pre- 
sented. As  we  have  previously  stated,  the  solu- 
tion of  this  problem  depends  only  on  the  value  of 
the  interaction  parameter  xt  and  the  initial  profile 
shape  factor  v. 

The  first  case  considered  is  for  xt  s 7. 5 and  i 
* 0. 5.  Figure  5 indicates  the  calculated  shock 
shape  and  sonic  line.  We  note  tha>  upstream  of 
the  interaction,  the  sonic  line  appe«  rs  near  the 
middle  of  the  boundary  layer.  In  addition,  in  the 
supersonic  flow  ahead  of  the  shock  and  outside  the 
boundary  layer,  the  flow  is  a simple  wave  type. 
Using  the  computed  values  of  <pj  along  the  line 
y */6*  = 1,  we  have  traced  the  characteristic 
waves  using  eq.  (7. 30).  It  is  seen  that  these 
compression  waves  emanating  from  the  curved 
sonic  line  within  the  boundary  layer,  intersect  the 
Impinging  normal  shock  causing  it  to  bend  for- 
ward. Although  the  characteristics  tend  to  merge 
together  near  the  shock,  there  Is  no  evidence  of 
shock  focusing  to  form  another  leg  of  the  shock 
pattern.  Also,  the  shock  blends  smoothly  into  the 
sonic  line  within  the  boundary  layer.  The  distri- 
butions of  the  normalized  perturbation  pressure  p 
defined  in  eq.  (7.31)  are  plotted  vs  S for  several 
values  of  y,  as  shown  in  Fig.  G.  This  figure 
gives  ar  indication  of  how  tlic  pressure  trans- 
forms from  a discontinuous  jump  far  from  the 
wall  to  a smooth  pressure  rise  at  the  surface. 
Comparison  of  the  distributions  at  y = 0,0.  5,  and 
1 illustrate  the  small,  but  significant  effect  of  the 
normal  pressure  gradient  through  the  boundary 
layer.  For  larger  values  of  y,  the  pressure 
steepens  to  a shock  wave  and  eventually  over- 
shoots the  undisturbed  value  at  y 4,  5.  For  guc- 
cocdingly  larger  values  of  y the  pressure  over- 
shoot decays  to  the  normal  shock  value.  The  gen- 
eral features  of  the  pressure  distribution,  includ- 
ing the  overshoots  arc  in  qualitative  agreement 
with  llu*  measurements  of  Gadd, ,,c'  for  a ease 
with  nearly  the  same  values  ot  x,  and  h.  How- 
ever, wo  note  tli.it  the  overshoots  in  Garid’s  data 
are  at  least  partially  due  to  an  axial  pressure 
gradient  induced  by  wall  effects  downstream  of  the 
shock  wave.  Another  feature  of  this  calculation 
is  shown  in  Ftg.  7.  Here  Ihc  normalized  velocity 
profile  u*  ul  - 1 •*  «(u,r  <>„)  with  e -0.034  is 
plolUd  vs  y for  several  values  of  x.  The  figure 
shows  the  velocity  profile  proceeding  from  the 
Coles’  law  of  the  wall-law  of  the  wake  form  at  up- 
stream infinity,  to  the  same  profile  uniformly  de- 
creased by  the  normal  shock  values  far  down- 
stream. The  intermediate  profiles  show  the  tran- 
sition between  the  two  extremes,  with  a smooth 
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transition  occurring  near  the  wall  and  an  abrupt 
change  far  from  the  wall.  The  velocity  at  x=  -0.  G 
•shows  the  profile  through  the  curved  shock?,  which 
it  intersects  at  approximately  y « 2. 5.  Altltough  it 
is  not  very  pronouuc  d in  the  figure,  the  velocity 
profiles  exhibit  over*  :»oots  down-stream  of  the  in- 
teraction at  several  xnmdary  thickness  above  the 
plate.  Overshoots  of  this  type  have  been  noticed 
in  experiments  on  shock  wave -bounce ry  layer  in- 
teractions (c.f.  Vidal  et  al,J7>). 


The  flow -field  structure  discussed  for  this 
first  case  is  typical  of  others  where  our  analysis 
is  appl  cable.  To  indicate  the  rffect  of  \t  on  our 
solution,  we  computed  results  at  y,-  5,7.  5,  10, 

20.  ,md  40,  all  with  r>  - 0.  5.  Figure  8 indicates 
the  effect  of  y,  on  the  normalized  surface  pertur- 
bati'in  pressure  distribution  vs  x.  We  see  that  in 
this  eooulinate  system,  the  pressure  rise 
steepens  for  increasing  y,  and  the  pressure  din* 
tiibuliui  becomes  only  slightly  more  asymmetrical 
w.’th  resjH-ct  to  the  incident  shock  position  (x  0). 

In  the  original  coordinate,  x-/};,  x,  the  scale  of 
the  interaction  is  magnified  and  the  asymmetry  of 
the  pressure  distribution  becomes  somewhat  more 
pronounced.  However,  we  do  not  sec  the  develop- 
ment of  a long  tail  down  .tream  of  the  shock,  as 
will  be  discussed  later. 


The  effect  of  5 on  the  initial  boundary-layer 
defect  profile  is  shown  in  Fig.  10.  Decreasing  » 
results  in  a more  “fuH*’  velocity  profile.  These 
profiles  were  then  used  to  compute  a set  of  outer- 
layer  solutions  at  Xt-  20.  From  Fig.  11,  we  see 
that  decreasing  $ to  achieve  fuller  initial  profiles 
results  in  steeper  surface  pressure  gradients. 


Next,  we  illustrate  our  numerical  solution  in 
comparison  with  the  free  interaction  solution  of 
Adamson  and  Feo<zu  (eq.  (B.  11)).  Here,  we  plot 
their  result  in  terms  of  our  variables  p and  x for 
the  case  yt  - 7.  5.  Their  solution,  however,  only 
gives  a surface  pressure  distribution  up  to  sonic, 
and  contains  an  arbitrary  additive  constant  to  tlic 
value  of  x.  Thus,  to  compare  it  to  our  solution, 
we  had  to  superimpose  Adamson  and  Feo’s  result 
and  attempt  to  match  the  pressures.  But,  com- 
paring to  our  surface  pressure,  the  agreement 
was  not  very  good.  However,  since  their  anal- 
ysts neglected  normal  pressure  gradients'through 
the  boundary  layer,  we  felt  it  to  be  more  appli- 


cable to  compare  to  our  results  at  the  edge  of  the 
boundary  layer,  y - 1.  The  comparison  shown  in 
Fig.  12  is  seen  to  be  quite  good  and  substantiates 
the  agreement  between  the  two  theories  for  sinal' 
Xf 


Next,  we  attempted  to  compare  our  calculation 
with  experimental  data.  We  had  to  restrict  our- 
selves to  normal  shock  data  from  a high  Kcynolds 
number  facility,  with  a fully  developed,  non-sep- 
a rated,  turbulent  boundary  layer,  '..*r  which 
=•  0(1).  The  experiments  which  came  nearest  to 
fitting  these  requirements  were  those  of  tiadd. 1361 
His  exjxTiment:*  included  a case  where  M„  1.  12, 
t‘  0. 0339  and  y,  7.502.  Unfortunately,  the  test 
took  place  in  a relative  .small  diameter  circular 
tube  with  strong  axial  pressure  gradients.  His 
post  shock  data  never  reached  the  n wise  id  shock 
jump,  even  at  largo  distances  away  from  the  wall. 
Nevertheless,  we  presojd  a comparison  of  our 
calculation  with  his  result  u Fig,  13.  Here  the 
notation  pVp®  corresponds  In  the  ratio  jf  static 
pressure  to  the  stagnation  pressure  upstream  of 
the  interaction.  Wo  see  that  the  numerical  and 
experimental  results  show  similar  overall 
trends. 


Fig.  9 indicates  the  variation  of  the  shock  wave 
and  sonic  line  as  a function  of  yt.  Wc  see  that  as 
the  interaction  parameter  increases,  the  sonic  line 
approaches  the  wall.  Also,  in  this  limit,  the  re- 
sults show  the  shock  wave  becoming  more  normal 
and  penetrating  further  into  the  boundary  layer. 


The  recent  experimental  data  of  Vidal  et  al, m> 
were  taken  in  a Ludwicg  tube  at  very  high  Rey- 
nolds number,  Re  36x10*.  Unfortunately,  their 
runs  were  at  M„  1.4  and  € - 0.  024  with  40.  5 
and  contained  a separation  bubble.  The  pressure 
seemed  to  have  an  upstream  influence  of  1-2 
boundary- layer  thickness  followed  by  a very  long 
tail,  (greater  Han  ten  boundary-layer  thickness) 
where  it  slowly  rises  towards  the  normal  shock 
value.  The  Mach  number  for  this  case  is  rela- 
tively high,  necessitating  the  use  of  the  extended 
small  disturbance  theory  (i.e. , eq.  (4.27)).  In 
addition,  the  data  shows  evidence  of  a weak  dis- 
turbance upstream  of  the  main  interaction.  The 
effect  of  this  disturbance  can  be  taken  as  a lower- 
ing of  the  interaction  Mach  number.  Hence,  for 
our  calculation  wc  used  a value  of  M„-  1.32  in- 
stead of  1.4,  for  which  the  similarity  parameter 
R , 7.5,  Again,  the  numerical  and  experimental 
results  show  reasonable  qualitative  agreement. 
?Ke  scale  of  the  upstream  influence  is  approxi- 
mately the  same,  but  the  downstream  scales  arc 
not  In  agreement.  We  note  that  this  calculation, 
using  the  extended  theory,  rcco  ered  the  exact 
pressure  jump  across  a normal  shock,  whereas 
the  experiment  showed  an  extremely  slow  ap- 
proach towards  the  normal  shock  value. 


Before  proceeding  with  a discussion  of  these 
results  in  the  next  l ection,  we  will  give  some  de- 
tails of  the  uumeri  al  calculations.  All  runs  con- 
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sidered  here  contained  78  mesh  points  in  the  4 
direction  and  57  in  the  tj  direction.  The  last 
mesh  point  before  infinity  in  the  physical  plane 
waa  at  x = ± 20  and  y = 20.  Typically  16  points 
were  included  within  the  boundary  layer  0 < y < 1 
with  a similar  mesh  in  the  inner  rectangle  (see 
Fig.  4)  in  the  x direction.  The  solutions  con- 
verge to  a rpaximum  residual  of  10**  in  100  - 800 
iterations  with  the  lower  value  corresponding  to 
X,=  40  and  the  higher  value  to  x»  = 5.  The  corre- 
sponding computation  times  were  from  1 to  6 min- 
utes on  an  IBM  370/168.  A relaxation  factor  u 
* 1.75  was  used  for  the  subsonic  region  and  u> 

« 1.0  for  the  supersonic  region.  A numerical  sta- 
bility problem  was  encountered  for  values  of  x,  . 

< 5,  for  which  we  could  not  obtain  any  converged 
results.  Another  difficulty  encountered  was  the 
use  of  the  wall  boundary  condition,  eq.  (7.  0). 

This  condition  was  apparently  unstable  and  we  had 
to  resort  to  the  original  wall  condition  eq.  (7. 6). 
We  circumvented  the  logarithmic  singularity  In 
My)  by  considering  the  wall  to  be  locate*',  at  y 
= 10**  instead  of  at  zero.  Numerical  studies  vary- 
ing this  small  wall  displacement  showed  negligible 
effects  on  the  calculated  surface  pressures.  Var- 
ious studies  also  were  made  on  the  effect  of  mesh 
size  and  mesh  distribution  on  the  accuracy  and 
stability  of  the  method.  To  illustrate  just  one  of 
these  results,  we  present  in  Fig.  15  two  calcula- 
tions for  X(-  10,  5=  0.  5 computed  with  8 and  16 
points  in  the  boundary  layer. 

DC.  Discussion 


Our  objective  has  been  to  develop  a rational 
approach  for  analyzing  the  interaction  of  shock 
waves  with  turbulent  Imumlary  layers.  Tin*  pres- 
ent analysis  considered  the  interaction  involving 
weak  normal  shock  waves  and  fully  developed 
turbulent  boundary  layers.  Under  those  condi- 
tions, the  sonic  line  is  located  in  the  main  part  of 
the  approaching  boundary  layer  and  as  a result  (ho 
suilace  pressur  rises  more  siu  /otlily  Ilian  in  in 
b taction.*:  at  higher  Ma'li  iminlr-rs.  We  have  «t<  - 
vc loped  a complete  asymptotic  theory  for  this  type 
of  in  I tract  ion  in  tin  coublt  limit  of  large  Rey- 
nolds numbers  and  Mach  iinmbeis  approaching 
one. 

The  asymptotic  structure  has  been  shown  to 
Involve  three  layers  and  we  have  determined  the 
leading  terms  of  the  solution  in  each  region.  The 
three-layer  structure  is  numerically  significant, 
since  the  middle  or  blending  layer  contributes 
logarithmic  terms  to  the  skin  friction  solution 
that  art*  multiplied  by  relatively  large  constants 
(l.c.t  r*/#c  *=  7. 3).  One  of  the  unexpected  findings 
of  the  present  study  was  the  discovery  that  the 


method  of  matched  asymptotic  expansions  fails 
for  compressible  turbulent  boundary  layers.  We 
have  shown  that  transformation  to  Van  Priest's 
generalized  velocities  circumvents  this  difficulty 
and  enables  us  to  obtain  complete  solutions  In  the 
inner  layers. 

The  asymptotic  theory  for  the  outer  region 
leads  to  a similarity  parameter,  Kt,  which  ia 
proportional  to  the  ratio  of  the  velocity  change 
across  the  shock  wave  to  the  velocity  defect  in  the 
upstream  boundary- layer  profile.  Examination  of 
the  outer  equations  indicates  that  the  solution  de- 
pends upon  the  value  of  two  parameters  Rt  and  i 
where  i is  a measure  of  the  profile  shape  in 
Coles'  correlations.  The  value  of  the  Reynolds 
number  and  Mach  number  influence  the  solution 
only  through  their  effect  onftjM.,  c}.  Thus,  a 
similarity  exists  for  moderate  strength  interac- 
tions (Kt  * 0(1)). 

We  have  demonstrated  that  our  analysis  for 
moderate  shock  strengths  reduces  to  Adamson 
and  Feo's,  au  theory  as  the  shock  strength  ap- 
proaches zero.  A numerical  comparison  has  in- 
dicated that  the  theories  are  in  good  agreement  ■ 
for  Xt  as  large  as  7.  5 (or  Rt  = 2. 125).  Our  results 
for  this  case,  given  in  Fig.  12,  show  a significant 
pressure  variation  through  the  boundary  layer, 
particularly  near  the  shock  wave.  This  pressure 
variation  Is  an  important  feature  of  the  flow  in  the 
region  whert*  the  shock  wave  penetrates  Into  the 
boundary  layer.  Nevertheless,  It  is  interesting 
that  Adamson  and*  Feo's  free  Interaction  solution 
is  apparently  quite  accurate  over  tlie  entire 
supersonic  region,  when  assumed  to  apply  at  the 
edge  of  tin*  boundary  layer,  rather  than  at  the 
wall. 


We  have  implemented  a mixed- flow  relaxation 
technique  for  solving  the  inviscid  outer-flow 
problem  and  have  presented  results  for  various 
values  t»l  Xi  between  5 and  40.  Numerical  con- 
vergence problem*:  prevented  us  from  obtaining 
solutions  for  values  of  x»  less  than  5,  Our  re- 
sults give  a reasonably  complete  picture  of  the 
structure  of  the  wave  pattern  in  a turbulent  shuck 
wave  interaction.  The  present  solution  and  the 
free- Interaction  solution  of  Ref.  21  indicate  that 
the  development  of  embedded  shock  waves  by 
focusing  of  compression  waves  In  the  upstream 
flow  does  not  occur.  Rather,  the  wave  pattern 
appears  to  develop  by  the  simple  mechanism  of 
compression  waves  being  generated  by  thickening 
of  the  subr.utlic  part  of  the  boundary  layer  and 
then  intersecting  the  approaching  shock  wave  and 
causing  it  to  bend  forward.  The  numerical  re- 
sults seem  to  indicate  that  the  shock  continues  to 


: 

I 


i 


21 


bend  forward  until  it  joins  smoothly  to  the  sonic 
line  in  the  boundary  layer.  Variation  in  the  pa- 
rameters Xt  and  i have  little  effect  on  the  general 
structure  of  the  wave  pattern.  However  the  value 
of  x<,  dot's  have  a strong  influence  on  the  scale 
of  the  wave  pattern.  As  increases  the  sonic 
line  approaches  the  wall,  and  the  shock  wave  be- 
ejincs  more  normal  with  most  of  the  interaction 
occurring  near  the  impingement  point  ol  the  un- 
disturbed shock  wave. 

Experiments  with  normal  shock  waves  .some- 
times indicate  the  presence  of  a supersonic  tongue 
and  a second  (bifurcated)  shock  wave  with  the 
overall  shock  structure  ri  sembling  a Mach  re- 
flection, A supersonic  tongue  is  evident  in  Sed- 
don,stM*  data,  but  not  in  the  data  of  Vidal  ert  al. 

The  present  numerical  results  do  not  show  the 
existence  of  either  a supersonic  longue  or  a sec- 
ond (bifurcated)  shock  wave.  The  reasons  for  the 
absence  of  these  features  in  Ihc  present  calcula- 
tion are  not  clear.  It  may  l>e  that  these  effects 
are  present  only  for  strong  shock  waves  involving 
separation  and  hence  arc  not  covered  by  the  sim- 
ple theory  developed  here.  It  is  also  possible 
that  the  present ‘numerical  scheme  is  not  accurate 
enough  »o  resolve  the  fine  details  of  the  How  near 
the  si  took  impingement  point.  A 'iner  mesh  and 
a more  accurate  finite-difference  method  than  the 
first-order  scheme,  used  for  supersonic  points  in 
the  present  calculation,  may  be  required  for  this 
purpose.  In  any  event,  the  present  results  indi- 
cate that  the  wave  pattern  is  not  a Mach  reflection 
in  the  usual  sense  of  the  term.  The  “forward 
limb.  ” described  in  experimental  results  is 
shown  to  be  a simple  continuation  of  the  main 
shock  wave.  The  "rear  limtf*,  when  it  is  present, 
is  just  a simple  recompression  of  the  flow  in  a 
supersonic  tongue  to  subsonic  speeds.  There 
does  not  seem  to  be  any  reason  for  assuming  that 
the  “rear  limb”  must  intersect  the  main  branch 
of  the  shock  wave  to  form  a triangular  shock  pat- 
tern. 

As  a final  remark  on  flow  structure,  we  note 
that  the  streamwise  length  scale  exhibits  a mini- 
mum as  a function  of  x».  The  minimum  length 
scale  occurs  for  x»  - 0(1)  and  increases  In  both 
limits  of  xt  large  and  small.  For  large  Xt  the  in- 
crease in  length  is  strongly  asymmetrical  with 
the  interaction  distance  increasing  on  the  down- 
stream side  of  the  shock  wave  and  becoming  van- 
ishingly small  on  the  upstream  side. 

The  prcaunt  theory  indicates  that  the  surface 
pressure  distribution  is  relatively  diffuse  and 
nearly  symmetrical  for  moderate  to  weak  shock 
strengths,  Xt < 10  and  becomes  more  one-sided 


and  steep  on  the  upstream  side  for  strong  shocks, 
Xt  > 20.  Although  this  general  trend  is  in  accord 
with  experimental  observations,  the  theoretical 
calculations  do  not  show  the  very  long  tail  char- 
acteristic of  surface  pressures  measured  behind 
normal  shock  waves  in  wind  tunnel  experi- 
ments. 

In  these  experiments,  the  pressure  distribu- 
tions downstream  of  the  shock  wave  showed  little 
indication  of  approaching  the  normal  shock  limit 
over  the  region  surveyed  by  the  experiment.  In 
3U  unit  37  the  pressure  was  measured  toalwut 
10  A*  behind  the  shock  wave  and  in  Ref.  30  the 
downsli  cam  extent  was  about  50  The  under- 
lying cause  of  these  long  pressure  “tails”  most 
likely  stems  from  wall  or  other  interference  ef- 
fects in  the  subsonic  flow  downstream  of  the 
shock  wave.  This  was  definitely  a factor  in  the 
measurements  of  Gadd.  °6*  He  called  attention  to 
the  presence  of  significant  axial  pressure  gradi- 
ents in  the  downstream  How  and  suggested  that 
they  were  due  to  displacement  effects  induced  by 
the  wall  boundary  layers.  The  experiments  of 
Seddon*18*  and  Vidal  ct  al<17>  contained  regions  of 
separated  flow.  The  presence  of  these  separation 
bubbles  could  have  a large  effect  on  the  down- 
stream How  and  be  a factor  in  the  formation  of 
pressure  tails.  Clearly,  further  experimental 
and  theoretical  efforts  are  necessary  to  under- 
stani  the  reasons  for  the  disparity  in  the  down- 
stream results.  However,  because  of  the  large 
regioi  of  subsonic  flow  behind  the  shock  wave.  It 
may  be  very  difficult  to  eliminate  or  correct  the 
data  for  wall  interference  in  experiments  on  nor- 
mal shock  waves  in  wind  tunnels. 

We  have  developed  a solution  for  the  leading 
terms  of  the  skin  friction  in  terms  of  the  outer- 
layer  wall  pressure  distribution.  Preliminary 
comparison  with  data  (not  presented  here)  Indi- 
cates that  this  solution  is  not  accurate  for  tran- 
sonic Mach  numbers.  The  difficulty  appears  to  be 
associated  with  the  anomalous  mean  dilatation  ef- 
fects, that  Bradshaw*26’  has  shown  to  be  impor- 
tant in  compressible  boundary  layers  with  pres- 
sure gradients. 

Although  the  skin  friction  is  not  accurate,  it 
ran  be  used  as  a guide  to  estimate  the  incipient 
separation  pressure  rise.  The  skin  friction  solu- 
tion ind  aus  that  the  change  in  skin  friction  is 
proportional  to  the  change  in  surface  pressure. 
Thus,  the  pressure  rise  required  to  separate  a 
turbulent  boundary  layer  is  0(1),  that  is,  indepen- 
dent of  Reynolds  number  in  the  limit  of  large  Rey- 
nolds number.  This  conclusion,  which  was  also 
arrived  at  by  Adamson  and  Feo,  CI*  is  in  broad 
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agreement  with  the  high  Reynolds  number  data  of 
Roskbo  and  Thomke. <m  It  is  interesting  that 
Townsend4*'  reached  this  conclusion  a decade  ago 
in  his  studies  of  turbulent  separation  in  low  speed 
flows. 

The  present  results  represent  the  first  step  in 
the  development  of  more  complete  theory  of  the 
interaction  of  shock  waves  with  turbulent  boundary 
layers.  The  results  suggest  the  following  prob- 
lems that  appear  worthwhile  for  future  study. 

• The  extension  of  the  present  theory  to 
strong  shock  waves  with  pressure  rises  of 
0<1).  Our  results  indicate  that  this  is  the 
range  of  shock  strength  that  is  important  for 
practical  problems  involving  separation. 

• The  development  of  higher-orde^  terms  in 
the  outer  expansions.  This  extension  should 
be  carried  with  a view  to  explaining  the  ab- 
sence of  the  long  pressure  tails  In  the  pres- 
ent theory.  The  understanding  of  the  pres- 
sure tails  is  essential  for  Improving  the 
agreement  between  theory  and  experiment. 

In  this  regard  it  would  be  very  useful  to 
carry  out  careful  experiments  for  interac- 
tions with  shock  waves  of  moderate  strength 
over  a range  of  x«  of  less  than  10.  High 
Reynolds  numbers  are  not  essential  for  this 
purpose.  However  special  care  should  be 
exercised  to  determine  the  effect  of  tunnel 
interference  on  the  development  of  the  pres- 
sure tails. 

• There  is  a need  to  Improve  llie  solution  for 
the  skin  friction.  The  present  results  indi- 
cate that  the  effects  of  anomalous  mean  dila- 
tation are  important  and  need  careful  con- 
sideration. 

• The  numerical  algorithm  lined  in  (lie  outer 
region  should  lie  improved  to  incluck*  a 
more  satisfying  implementation  of  the  wall 
boundary  condition  and  lo  obtain  greater  ac- 
curacy ir  the  supersonic  region.  There  is 
a need  to  develop  a second-order  accurate 
Unite -difference  scheme  lor  .*n  per  some 
point:  in  a mixed  flow  relax. tie  i:  procedure. 
This  woJrl  allow  greater  loolmim  near 
the  shock  impingement  point  and  enable  a 
more  careful  study  ot  "Mijh-rsoiiie  tongues*' 
and  tlK*  second  shock  wave.  U would  also 
be  useful  to  uncover  the  reasons  for  (ire 
lark  of  convergence  of  the  iterative  solution 
f »**  small  Xt.  Wi  < 5). 

• The  present  theory  is  applicable  only  to  flat 
surfaces.  T he  effects  of  even  a small 
amount  of  cm  ■•ure  arc  known  to  be  signif- 
icant. Inviscid  : eory  indicates  the  pres- 
ence of  a sirring  . ugular  expansion  just 
downstream  of  a norma]  shock  wave  on  a 


curved  surface.  Experiments  on  airfoils 
indicate  that  the  pressure  behind  the  shock 
wave  differs  from  the  normal  shock  pres-  . 
sure  determined  from  the  shock  j*«mp  con- 
ditions. The  data  shows  that  pressure  be- 
hind the  shock  wave  Is  usually  quite  close 
to  sonic,  independent  of  shock  strength. 

The  reasons  for  this  are  unknown,  although 
it  is  clear  that  it  must  be  related  to  a com- 
plex interplay  between  the  inviscid  expan- 
sion and  the  shock  wave-boundary  layer  in- 
teraction process.  Thus  it  is  clear  that  the 
present  theory  must  be  extended  to  include 
surface  curvature  effects  if  it  is  to  apply,  to 
interactions  of  practical  importance.  , ’* 
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Appendix  A,  Matching  of  the  Non-interacting 
Velocity  Profiles 

In  this  appendix,  we  consider  the  Maise-Mc- 
Donald*80’  profile  given  by  eqs.  (3. 1)  through  (3. 5) 
as  an  exact  solution,  and  show  that  the  usual  limit - 
functicn  expansions  do  not  overlap.  Expanding 
(3. 1)  and  (3. 4)  in  powers  of  c with  y fixed,  we  ob- 
tain 

“•«urcl  + «ui-€t2*nX  + ***.  (A-1) 

where 

«i=“logy -~(l+coswy)  (A. 2) 

The  inner  (limit-function)  expansion  is  obtained  by 
expanding  (3. 1)  and  (3. 5)  for  small  € hold’ng  y* 
fixed,  whereby 

(a.s) 

with 

u*,  = ^ logy*  + B0  (A.  4) 

In  order  to  demonstrate  matching,  we  must 
compare  the  inner  expansion  of  the  outer  solution 
with  the  outer  expansion  of  the  inner  solution. 
Carrying  out  these  constructs  and  expressing  both 
expansions  in  the  outer  variable  y,  we  obtain 

Inner  [uowtvr1:  1 4gl  + 0(c*)  (A. 5) 

Outer  fuIM„]  (u,  - aii*  4 • • • ) + c (1  - iawj 

♦ •••)g,  + 0(e*)  (A. 6) 

where 

g^-logy--  (A.  7) 

The  coefficients  of  cq.  (A.G)  can  ho  shown  to  bo 
an  infinite  series  in  power::  of  (au,)  arising  from 
contributions  of  higher  order  terms  in  the  Inner 
solution.  Comparison  of  eqs.  (A. 5)  and  (A.G) 
shows  that  the  inner  and  outer  limit  function  ex- 
pansions do  not  match  for  non- zero  M„.  Analysis 
using  the  more  general  notation  of  intermediate 
variables,  (o.g.  sec  Kef.  32)  indicates  that  this 
difficulty  cannot  bo  rectified  since  the  domain  of 
validity  of  the  two  expansions  does  not  overlap. 

The  matching  problem  can  be  viewed  as  being 
caused  by  the  simple  expansion  of  the  trigonomet- 


ric terms  in  the  inner  solution.  Such  an  expan- 
sion implies  that  the  density  is  a constant,  equal 
to  the  wall  value,  in  the  lowest-order  solution. 

The  problem  Is  resolved  by  assuming  that  the  so- 
lution in  the  wall  It  yer  Is  given  by  the  full  expres- 
sions In  eqs.  (3. 1)  and  (3. 5).  We  now  show  that 
the  generalized  inner  expansion  matches  the  outer 
limit-function  expansion  given  by  eq.  (A.  1). 

To  demonstrate  matching,  we  write  the  full 
wall  layer  solution  In  outer  variables 

«i—r=|sina|u.4^=[^logy (l-cos»y)]j 

(A.  8) 

We  assume  this  solution  is  valid  In  a domain  y* 

- «>  such  that 

0(c  logy  )<  0(1) 

This  allows  us  to  expand  Eq.  (A.  8)  for  small  (logy 
which  results  in  the  following  large  y*  solution:  t 

Outer[utiMr]  = 1 +(Uj-  (*~,uf  (A.  9) 

Comparison  of  this  last  result  with  eq.  (A.  1) 
demonstrates  that  the  generalized  inner  solution 
matches  the  outer  (limit-function)  expansion. 

This  demonstrates  that  the  Maise-McDonald 
correlation  does  provide  an  initial  profile  that  is 
consistent  with  the  asymptotic  formulation  of 
compressible  turbulent  boundary  layers.  In  this 
analysis,  we  assume  the  initial  profile  to  be  given 
by  eqs.  (3. 1)  through  (3.  5).  This  will  allow  a ra- 
tional iwo  parameter  fit  of  specified  experimental 
profiles.  We  also  note  that  the  profiles  can  in* 
represented  by  asymptotic  expansions  given  by  Hq. 
(A.  1)  in  the  outer  region  and  by  eq.  (3. 5)  with  the 
5 terms  neglected  in  the  fully  turbulent  part  of  the 
Inner  region. 

Appendix  It.  Solution  of  t lie  Outer  Layer  for  0 

The  outer-layer  equations  given  by  eq.  (4. 14) 
should  reduce  to  Adamson  and  Feo's<ll>  weak 
shock  theory  in  the  limit  Xt**0.  The  approach  to 
this  limit  is  somewhat  complicated  and  involves 
the  splitting  of  the  outer  inviscid  region  into  two 
inviscid  layers.  One  layer  consists  of  the  main 
part  of  the  defect  region  and  is  scaled  by  the 
boundary-layer  thickness.  The  strcamwisc  scale 
must  be  stretched  by  x,  as  follows 

(B.l) 
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The  y coordinate  remains  unchanged.  Thus,  inner 
variables  are  given  by  il,  and  y.  Eq.  (B.  1 ) indi- 
cates that  the  streamwise  length  scale  increases 
as  xt  gets  small  and  this  implies  that  normal 
pressure  gradients  in  the  defect  layer  will  vanish 
to  lowest  order  as  Xi~0.  In  this  inner  region,  <t> 
should  be  expanded  In  the  form 

(u.2) 

Substituting  eqs.  (B.  1)  and  (B.  2)  into  eq.  (4. 14) 
and  imposing  the  boundary  condition 

<MM)=o 

leads  to  the  following  solution  #lr: 

yto  id- 3/ 


This  is  the  transonic  small  disturbance  equation 
which  Adamson  and  Feo  showed  to  govern  their 
outer  solution.  Equation  (B.  9)  must  be  solved 
subject  to  the  boundary  conditions 

£ t - x(y )dy  = 0 for  $ - 0 

(B.  10) 

<?i-0forx 

which  is  derived  by  matching  the  Inner  and  outer 
solutions.  The  problem  formulation  is  completed 
by  specifying  boundary  conditions  for  large  y. 

The  boundary  value  problem  consisting  of  eqs. 

(B.  9)  and  (B.  10)  and  the  outer  boundary  conditions 
determined  by  the  incoming  wave,  is  precisely  the 
problem  arrived  at  in  Ref.  21.  They  showed  a 
turbulent  free  Interaction  solution  can  be  derived 
by  using  a simple  wave  solution  of  eq.  (B.  9). 


where  at  this  stage  of  the  solution  is  an  arbi- 
trary function  and  X is  defined  in  eq.  (4. 15).  The 
above  solution  Implies  that  the  pressure  is  con- 
stant in  the  defect  layer  and  is  given  by 


fc-MO—  *o»  (B-4> 

Thus,  eq.  (B.  3)  provides  a relationship  »v»tween 
the  normal  component  of  velocity  and  th-  pressure 
gradient 

(B'5> 

cq.,(B.  5)  is  the  defect  layer  solution  derived  in 
Ref.  21.  Notice  that  Vt  approaches  a non- zero 
limit  as  y approaches  Infinity, 

V,(*. «)  = V,.  = - x(y)dy  (B.  6) 


In  the  notation  of  the  present  paper,  their  free 
Interaction  solution  can  be  written  in  the  form 


*.*-*-,[* -M 

(B.  11) 

•A“"CJa-)l 

where  A is  the  constant  appearing  in  Ref.  21  and 
is  given  by 


A * - § />+!  f*  x{y)dy 
or  using  eqs.  (3. 13)  and  (4. 15) 


(B.  12) 


(B.  13) 


Thus  the  Inner  solution  cannot  satisfy  the  outer 
boundary  conditions  and  an  outer  solution  is  re- 
quired for  large  y.  Analysis  indicated  that  the  ap- 
propriate outer  variable  is 


The  free  Interaction  solutions  given  In  eq.  (B.  11) 
form  a one-parameter  family  of  solutions  corre- 
sponding to  an  arbitrary  shift  of  the  origin. 


y=*/x»  (b.7) 

In  order  to  match  the  defect  solution  in  eq.  (B.  6), 
we  must  stretch  the  potential 

<P  = ^ 0 (B.  8) 

Substitution  of  outer  variable  ft,  y,  C>  into  eqs. 

(4. 14)  and  setting  Xt t0  zor°  r*»s*ult«  in  the  follow- 
ing equation: 

U+(r»D#,His-<i»*o  (b.9) 


The  solution  given  by  eq.  (B.  11)  leads  to  a 
surface  pressure  distribution  that  rises  smoothly 
from  free  stream  pressure  at  negative  infinity  to 
sonic  pressure  at  the  downstream  terminus  of  the 
solution.  The  pressure  gradient  rises  monotoni- 
cally  from  zero  at  negative  Infinity  to  a finite  val- 
ue at  the  sonic  point.  The  second  derivative  is 
positive  over  the  entire  interval  of  definition  of 
the  solution.  From  this  last  property,  it  can  be 
shown  that  the  simple  wave  pattern  does  not  in- 
volve a focusing  of  characteristics  nor  a forma- 
tion of  an  imbedded  shock  wave.  It  also  follows 
that  the  wave  pattern  is  terminated  downstream  by 
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One  of  the  principal  consequences  of  the  large 
Reynolds  number  limit  Is  the  fact  that  the  total  1 

shear  stress,  laminar  plus  turbulent,  is  constant  , 

to  all  orders  in  c across  the  wall  layer.  Analysis 
of  the  error  terms  in  the  wall-layer  equations  in* 
dicate  that  the  main  error  terms  for  large  y* 
arise  from  variations  of  total  stress.  It  follows 
from  consideration  of  these  terms  <hat  the  gen- 
eralized wall-layer  solution  is  valid  in  a region 
defined  by 

c-0  ) 

> such  that  0fcH(Y)]<  0(1)  (C.5) 

y*— ) 

It  follows  from  this  result  that  the  function  G 
defined  in  eq.  (C.4)  vanishes  for  c-0  in  this  re-  , 

gion.  Hence,  an  outer  expansion  of  the  inner  so- 
lution, valid  in  the  region  defined  by  eq.  (C.  5), 
can  be  obtained  by  expanding  eq.  (C.  3)  for  small 
G.  This  yields,  to  terms  of  second  order 


an  envelope  (or  caustic)  of  the  straight-line  char- 
acteristics. The  envelope,  which  is  convex  to  the 
supersonic  region,  starts  off  normal  to  the  x-axts 
at  the  sonic  point  and  is  asymptotic  to  a straight  line 
inclined  at  the  free-stream  characteristic  direc- 
tion at  positive  infinity.  The  implications  of  these 
properties  of  the  frec-interaction  solution  on  the 
wave  pattern  in  shock  wave-boundary  layer  inter;* 
actions  are  discussed  in  section  IX  of  this  paper. 

Appendix  C.  Matching  of  the  Blending  and 
Wall-Layer  Soiutions 

The  wall-layer  solution  using  Van  Driest's 
generalized  velocities  will  be  shown  to  match  the 
solution  in  the  blending  layer.  The  requirement 
that  the  two  expansions  match,  leads  to  a corre- 
sponding solution  for  the  skin  friction,  which  will 
be  derived  in  this  appendix.  Since  the  wall -layer 
solution  Is  not  in  the  form  of  a limit-function  ex- 
pansion, matching  of  the  two  expansions  requires 
some  special  steps.  First,  we  write  the  general- 
ized wall-layer  solution  given  by  eq.  (6. 27)  in 
terms  of  the  blending- layer  variable  y.  The  def- 
inition of  the  stretched  variables  y and  y*  and  the 
skin-friction  law  for  the  non-interacting  flow  can 
be  used  to  establish  the  following  relationship  be- 
tween y and  y*: 

ilo8y,.D0  = ’/Ssii.+  |;lOBC.H(Y)  (C.  1) 

where  I!(Y)is  the  funclion  defined  in  eq.  (5. 14). 
Equation  C.  1 is  used  to  express  all  y*  terms  ap- 
pearing in  the  wall-layer  solution  for  u,  eq. 

(G.  21)  intei  ms  of  the  blending-layer  variable  Y. 
The  result  can  be  written  in  ihe  form 

u - - sinjau,  i jp-  (J(Y;« ) J (C.  2) 

or,  expanding  the  trigonometric  terms  and  us  in;; 
tic  didinith'M  of  u, 


(C,  3) 

The  function  G is  given  by 


Outer  (Blending- Layer)  (uUw]  = 1 + |^clog< 

♦ <■'*“;«  T.!i  ♦ £(r.i  - P.i) 

‘ 2(r"' " 

+ i(m,q-l)  (rwl-pwl)IKY)+ 3(i  + m#q) 

* frwi  - Pwi)8  + a B0(r - pwl) 

“ 2^tw*  " P*r4  Pwi>j  (C.  6) 


The  previous  result  is  to  be  compared  with  the 
inner  expansion  of  the  blending-layer  solution. 
This  is  obtained  by  substituting  eq.  (5.28)  into 
eq.  (5.  7)  which  yields: 

3 

Inner  [u  l + y-*  log* 

. ■ (H(Y)-  P.l-eios '^(U 

. r’logc  [<1 0, n.)(|;)  P.  - ^7  H(Y)] 

- • c'[y HHY)-  J(t„ - M'.P,)H(Y)t  P, - C2] 

(C.7) 

Analysis  of  error  terms  indicate  that  this  expan- 
sion is  valid  in  the  domain  defined  by 
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C ~clog€  i c {l!(Y)+  i(r w| - pwl)q} 

CK 

+ |z  1 ~ (rvl  - pwl)| 

♦ c2|i  frwj  - Pwj)  (11  - »o) 4 2 

+ P*f-i(rwl-pw,)2]J  (C.4) 
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such  that  0(<  log  Y)  < 0(1) 


This  can  be  shown  to  overlap  with  the  region  of 
validity  of  the  wall  layer  defined  by  eq.  (C.  5). 
Thus,  eqs.  (C.  6)  and  (C.  7)  can  be  matched  term 
by  term  and  the  introduction  of  intermediate 
variables0*’  is  not  required  for  this  problem. 
Taking  note  that  the  pressures  match,  the  re* 
quirement  that  eq.  (C.  6)  and  (C.  7)  match  yields 
the  following  expressions  for  the  first  three  terms 
of  the  skin  friction  solution: 


r"4Ml-q)P 


(C.  9) 


Substitution  of  these  results  in  eq.  (6. 22)  leads  to 
the  skin-friction  solution  given  in  eq.  (G.  28). 


We  note  that  the  expansion  for  the  skin  fric- 
tion, eq.  (6. 28),  contains  a sequence  of  terms 
that  proceed  in  powers  of  the  function  (eyMip^j). 
Because  of  these  terms,  the  formal  expansion  for 
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the  skin  friction  approaches  its  large  Reynolds 
number  limit  very  slowly,  and  Is  useful  only  for 
very  small  free  stream  Mach  numbers  (e.g.  M* 

< 0. 5).  These  terms  arise  from  the  expansion  In 
powers  of  c of  the  wall  density  written  as  a func- 
tion of  the  surface  pressure.  Thus,  they  arise 
only  in  compress ibie  flows  with  pressure  gradi- 
ents. The  utility  of  the  formal  solution  for  the 
skin  friction  can  be  improved  by  a Euler  summa- 
tion of  these  terms.  The  rearranged  solution  can 
then  be  written  in  the  form 


(C.10) 


(i  + CyM Up., i cpwt + . . . ] K 


In  Uk  above  rearranged  series,  the  wall  density 
terms  are  effectively  summed  to  all  orders  in  c. 
Notice  that  contributions  from  the  wall  density,  in 
the  second  term  of  the  above  solution,  now  vary 
between  one  for  M . = C and  zero  for  M.-®.  Pre- 
viously, these  terms  b came  unbounded  as  the 
Mach  number  approacm  1 infinity. 
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